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Abstract 

The content of arsenic, antimony, selenium and tellurium in soils from a territory with 
natural and man-caused load, located in the range of the Old Quaternary proluvial plume in 
the northeastern part of the Sofia valley, was studied. The contents have been determined by 
LA-ICP-MS, together with those of more than 45 chemical elements, for the purpose of 
future remediation and study of the impact of natural and man-made factors on the 
geochemical relationships of the elements in the specific soil environment. The studied soil 
types are Regosols and Cambisols, most of which have been changed and classified as 
Technosols (USS Working Group WRB, 2015, Teoharov et al., 2019, Stoykova, 2021) under 
the influence of industrial waste from the former Kremikovtzi metallurgical plant. Very high 
contents of the groups of elements Te - Se and As - Sb have been established, which are 
inherited from the soil - forming rocks with predominant, respectively syenite or ore material. 
The content of arsenic in Technosols is on average 4 times higher than the maximum 
permissible concentration of metalloid in industrial soils, and the arsenic in the control profile 
of Cambisols exceeds 1.5 times the norm for the maximum permissible concentration for 
arable land. The processes of overwetting, acidification and microbiological activity are the 
factors that influence the redistribution of the four elements in the soil environment. 
Contamination with agglomeration dust leads to the accumulation of antimony in the spolic 
horizon, where it is associated with a number of heavy metals. Its more metallic character is 
reflected in the statistically significant high positive correlation with iron - r (Sb-Fe) = 0.85, 
which shows only this element. 

Keywords: arsenic, antimony, contamination, Cambisols, Regosols, selenium, 
tellurium, Technosols 


Introduction 

Arsenic (As), antimony (Sb), selenium (Se) and tellurium (Te) are chalcophilic 
elements of the p-block of the Periodic Table, satellites (As, Sb) or analogues (Se, Te) of 
sulfur (Stefanova, 2005). In the primary endogenous environment they are concentrated in 
sulfides, related minerals and sulfosols, and in the exogenous environment, like sulfur, they 
are oxidized relatively easily (Stefanova, 2005, Plant et al. 2003, Fersman, 1939) and in the 
form of variously charged oxygen radicals compounds are involved in active soil metabolism. 
Soil macro and microbiome, including man, are directly influenced by the excess content of 
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these elements (Sadiq, 1995, Denys, 2009, Missen, 2020, Cuvardic, 2003). Their study is of 
interest to soil science not only because of their toxic effects on living organisms, but also 
because of the food needs of selenium, which is a trace element essential for organisms. 

In the context of the specific geological and ecological situation on the territory of the 
former Kremikovtzi MP and beyond, the study of these four elements serves purely scientific 
interests, as processes can be identified - natural and technogenic, affecting their distribution 
in the soil environment. The soils subject to the present study were developed on proluvial 
rocks containing materials of diverse geological origin, naturally enriched in metals and 
metalloids, primarily gold (Stoykova, 2021). The proluvial rocks, shaped like proluvial cones 
connected in an extensive proluvial plume, are the result of drainage through ravines of two 
feeding provinces (tentatively called western and eastern) that operated during the 
Pleistocene. They are poorly fused gravels with sandy-clay matrix and alternating grains of 
sand, gravel and clay horizons. They contain rock fragments and deposited products of 
Paleozoic argillites, siltstones, lydites and others. with varying degrees of metamorphism, 
carbonate and clay-carbonate Triassic and Jurassic sediments (Angelov et al., 2008), 
potassium-alkaline magmatic rocks from the Buhovo-Seslavtzi pluton (Dyulgerov, 2005), ore 
rocks and ore minerals from the Buhovo uranium deposit (Kalaydzhiev, 1993). Due to their 
loose structure and proximity to the mountain slope, they are relatively easy to weather and 
should have lost much of their primary minerals. According to Stoykova (2021), the 
technogenic situation in the district and the reported alkalization also invariably influenced 
the geochemical relations of the elements sensitive to the parameters of the environment, 
such as arsenic, antimony, selenium and tellurium. 

The complicated technogenic situation is due to the nearly 50 years of operation of the 
former Kremikovtzi metallurgical plant, which is a conglomeration of heavy industry plants, 
located on an area of 24,000 m?. Agricultural land around the plant has been the subject of 
many studies over the years. From them it can be concluded that the lands have high levels of 
heavy metals and arsenic around the plant, and with the removal of it only lead and arsenic 
remain exceeded (Faitondjiev et al., 2000; Shulin et al., 2007). 

This publication is the only one that interprets the results for the content of arsenic, 
antimony, selenium and tellurium, as a group obtained together with data for more than 45 
other chemical elements in soil profiles located near various factories inside the plant. 
External control profiles were also included in the study. 


Materials and Methods 

Soils from 15 profiles in the area of Kremikovtzi MC have been studied, incl. 11 from 
the territory of the plant - with man-caused influence and 4 controls in different directions 
outside it. The profiles are divided into a total of 6 geographical positions. In turn, the eleven 
profiles with man-made impact are grouped into three groups depending on the type of 
pollution, as well as the local geographical location: 

o In the middle parts of the western proluvial cone, in a border area with modern 
proluvium, rich in deposited limestones and potassium-alkaline magmatic rocks, mainly 
syenites (Terrain group I) - Profiles 1, 2,3,4+ 13c, 14c, 

o In the upper parts of the eastern proluvial cone, rich in deposited limestone and ore 
matter - (Terrain group II) - Profiles 5, 6, 7, 8. 
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o In the middle parts of the eastern proluvial cone - (Terrain group HI) - Profiles Ne 9, 
10, 11. 

A number of standard soil analyzes were performed on the 15 soil profiles. In view of 
the large database, five man-made influences were selected (Profiles Ne 3, 6, 7, 8 and 11) and 
one control profile - Ne 13c, the most representative for each of the initially separated groups, 
which were subjected to additional chemical and mineralogical research, incl. LA-ICP-MS. 
The results of the chemical analyzes are an interpretation based on the parameters 

o Average content, 

o Linear Brave correlation coefficient (r). The statistical significance of the obtained 
coefficient is calculated by the formula | 7 | > 3 * or, where or = (I - r°) /Nn, and “n” is the 
number of pairs in the correlation. After the calculations, it was determined that the 
correlation coefficients in the intervals (-1 + -0.5) and (0.5 ~+ 1) (Rollinson, 1993) were 
statistically significant. 

o Concentration coefficient (CC), calculated as the ratio between the content of the 
element in the specific horizon of the profile (Ke) to its content in the lowest horizon of the 
same profile (Kc) — Ke / Kc (Fersman, 1953; Stefanova, 2005). This coefficient makes it 
possible to identify accumulation and migration processes vertically, within the soil profile. 

According to their genetic nature and the natural course of soil development, 
secondarily contaminated soils are divided into two groups: Cambisols - profile 3, profile 7, 
profile 13c and Regosols - profile 6, profile 8, profile 11. According to the changes that have 
taken place due to the influence of industrial materials, some of them are already classified as 
Technosols (Stoykova, 2021). 

Morphologically, the studied soils are reddish in color (50% of all samples have a hue 
of 5 YR), with high fragmentation, naturally or artificially overlapped. Cambisols have well- 
formed genetic horizons, and Regosols have a layered profile structure. Soil-forming 
sediments are reached everywhere - alternating gravel, sand and clay layers. Vegetation in the 
technogenic profiles is scarce or absent, except for one secondary afforested with acacia 
profile (Ne 1). 

Results and Discussion 

According to the results obtained, the average content of the discussed elements in the 
studied soils is much higher, compared to selected reference values (Kabata-Pendias, 2011) - 
Te (56x) >> As (25x)> Se (24x) >> Sb (10x). The determined average contents are (in pg / g) 
As — 170,66 (409.34 + 24.63); Se — 10.56 (14.25 + 6.15); Sb — 6.66 (13.93 + 3.13) and Te — 
4.49 (7.03 + 2.14) (Table 1). 

Of the four elements, the content of arsenic in soils is subject to state regulation, on 
the one hand as the most common of them, but also because of its anthropogenic contribution 
to the soil and its toxic properties, falling into it. The values obtained for this element exceed 
the trigger values from Regulation 3 / 1.08.2008 as follows: arsenic in Technosols is on 
average about 4 times above the maximum permissible concentration of metalloid in 
industrial soils (40 ug / g, Regulation 3 / 1.08.2008), and the arsenic in the control profile 
exceeds 1,5 times the established norm for the maximum permissible concentration in arable 
land (25 ug / g, Regulation 3 / 1.08.2008). 
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Table 1. Average contents of the studied elements, compared to world reference 
values (Kabata-Pendias, 2011), distributed by geographical area and by man-made 
pollution. 


Average | | Avarage Field | Field | Field ey 
in soils,| |in  this|Mediana |Max | Min group |group | group a ne 
ref. data | | study I II Il ; 
avarage 
Lg/g 
As |6.83 170.66 |150.22 | 409.34 | 24.63 97.80 | 247.56 | 144.93 ee APA | 
Se |0.44 10.56 11.09 14.25 |6.15 11.82 |10.75 | 7.26 11.56 
Sb | 0.67 6.66 6.23 13.93 |3.13 8.27 |7.02 |5.89 4.04 
Te |0.08 4.49 4.95 7.03 |2.14 5.58 |4.43  |2.91 4.81 


The content of Te and Se is inherited from proluvial rocks enriched in syenite 
fragments, and that of As and Sb is bound to soil-forming materials enriched in ore minerals 
and sulfosols, arsenides and antimonides, some of the minerals accompanying uranium 
mining (Kalaydzhiev, 1993). This statement is based on the higher measured values for the 
two groups of elements in the profiles of the respective terrain groups, which are formed on 
different soil-forming rocks - Te and Se in profiles NeNe 3 and 13c of terrain group I and As 
and Sb in profiles NeNe 6, 7, 8, 11, part of terrain groups II and HI (Table 1, Table 2). The 
mechanism of Te concentration in alkaline magmatic rocks is known to science. Such rocks 
are associated with the largest deposits of tellurium and gold, in the form of tellurides, which 
are generally part of the ore (sulfide) fraction of accessory minerals (Kostov, 1993). 
Selenium, following its chalcophilic nature, is also concentrated in this fraction, where it 
forms its own selenides or replaces sulfur in metal sulfides (Kostov, 1993). 

Natural and technogenic soil-forming processes are imposed on these soil-forming 
rocks, which lead to redistribution of the studied elements within the soil profiles. 

The over wetting process covers the entire volume of profile 8 (Reductic Technosols, 
colluvic), combined with a strongly acidic soil reaction in the intermediate horizons (Table 
2). 

This has mobilized all four elements (CCTe = 2.41, CCSb = 1.55, CCSe = 1.86, 
CCAs = 2.43) (Fig. 1), which are known to be reduced in a reducing environment. easily 
reduce to neutral or weakly negatively charged mobile forms such as HAsO)”, H»AsO*, 
HSbO,”, SbO”, Se (s), Te (s) (Sadig, 1995, Denys, 2009, Missen, 2020, Plant et al., 2003) 
and hydrogenselenide anions (HSe’) - and hydrogenselenide (H2Se°), as a manifestation of 
the halogen nature of selenium (Cuvardic, 2003, Plant et al., 2003). 

Antimony is an exception to the general trend in this example. The upward 
movement is to a lesser extent, and fixation is observed in the most acidic horizon ACyy (125 
- 225 cm) and CC (1.55 and 1.53) for the surface and intermediate horizon, respectively. The 
behavior of antimony under these conditions reflects its metallic nature (expressed in the 
lowest energy of positive ionization Z“ among the four elements), by virtue of which in 
oxidation conditions it does not form charged oxygen radicals but solid mineral phases - 
Sb20,4, Sb2Os ( Stefanova, 2005). They are adsorbed and incorporated isomorphically into the 
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poorly soluble Fe-Mn oxide and hydroxide structures crystallizing during the drought and 
oxidation period, thus gradually accumulating (Denys, 2009). 

Table 2. Content of the studied elements, iron, manganese, sulfur and values of the 
main soil indicators. 


Base-cation 
saturation, | Total Physical 
Mn CEC % Carbon, | Carbonates, | clay, 
mceg/g cmol/kg | % CEC % % % 
Profile 3 Ekranic Spolic Technosols (Alcalic, Cambic, Chromic, Endoskeletic, Eutric, Toxic) 


Horizon / Depth, cm 


Asp 30-0 very black mass, agglomeration powder 


Aek 25-0 concrete flooring 

Aspok 0-50 13.25 | 9447.42 | 2726.75 | 79.10 | 10.35 | 13.93 | 4.68 
50-105 4.17 | 2117.42 | 1824.95 | 137.53 | 12.70 | 6.14 5.73 
105-170 
170-235 | 3.00 | 1239.74 | 1202.65 | 63.56 | 10.31 6.43 5.38 
235-347 | 4.65 | 1859.27 | 1872.34 | 111.01 | 13.92 | 6.58 6.52 


Profile 6 Spolic Technosols (Alcalic, Chromic, Colluvic, Epicalcic, Eutric, Toxic) 


Asp, 30-0 granular and non-granular coke powder 

Aspok 0-40 4.28 | 9333.53 | 2383.79 | 170.70 | 11.40 | 6.31 4.83 
40-60 5.94 | 7918.97 | 1402.77 | 322.38 | 11.23 9.22 5.76 
60-110 
110-160 | 5.21 | 4014.86 | 1534.74 | 275.37 | 10.94 | 7.78 4.87 
Cin g 160-260 | 6.35 | 13075.71 | 2393.49 | 409.34 | 11.42 | 8.93 5.09 


Profile 7 Spolic Technosols (Cambic, Chromic, Eutric, Gleyic, Toxic ) 
Asp 10-0 coal dust 
0-38 4.73 | 6042.74 | 2515.69 | 196.01 12.58 
38-68 
68-110 4.96 | 2736.36 | 1112.65 | 162.90 | 8.85 
110-165 
165-210 4.58 | 2737.84 | 1695.98 | 227.42 12.53 


Profile 8 Spolic Reductic Technosols (Chromic, Colluvic, Eutric, Gleyic, Toxic) 


Asp 20-0 tar. tar fus. naphthalene 

AC; g 0-40 5.36 | 3409.46 | 1728.74 | 333.90 | 13.34 | 7.85 5.51 
ACy 40-125 
ACin 125-225 | 5.09 | 1841.76 | 941.62 | 240.12 | 8.02 2.16 
ACw g 225-400 | 4.91 | 2189.57 | 957.05 | 137.42 | 7.17 2.29 


Profile 11 Spolic Technosols (Chromi,. Colluvi,. Eutric, Skeletic) 
Asp 20-0 loose sandy quartz substance 


Asp2 20-0 dark cast iron - slag 

Asp3k 0-30 5.26 | 4134.00 | 1443.99 | 108.41 
30-115 4.89 | 2489.79 | 1178.28 | 126.05 
Cu 115-215 5.82 | 2854.41 987.06 | 200.32 


Profile 13c Rhodic Eutric Cambisols (Clayic) 
0-32 3.72 894.48 811.88 
32-67 
67-116 3.88 1037.50 | 1502.27 
116-151 
151-185 
185-220 ‘ 1652.07 
220-260 
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Table 3. Concentration coefficients (CC) of the investigated elements, calculated as 
the ratio between the values in each of the horizons to the values in the lowest horizon of the 


profile. 
CC Te Sb Se AS 
Prof. 3 
Aspok 0.72 2.12 0.74 0.71 
B, 0.88 0.93 0.91 1.24 
1C 0.83 0.98 0.74 0.57 
2Cg 1.00 1.00 1.00 1.00 
Prof. 6 
Aspok 0.95 0.71 1.00 0.42 
aC; 1.13 1.03 0.98 0.79 
Cr 0.96 0.87 0.96 0.67 
Cin g 1.00 1.00 1.00 1.00 
Prof. 7 
Ak 1.27 1.00 1.00 0.86 
Bw(g) 0.42 1.10 0.71 0.72 
C 1.00 1.00 1.00 1.00 
Prof. 8 
AC, g 2.41 1.55 1.86 2.43 
AC 0.94 1.53 1.12 1.75 
ACy g 1.00 1.00 1.00 1.00 
Prof. 11 
Asp3k 1.32 0.61 1.36 0.54 
AC; 1.11 0.58 1.19 0.63 
Cr 1.00 1.00 1.00 1.00 
Prof. 13c 
A 0.34 0.66 0.49 0.49 
Bw 0.72 0.58 0.95 0.72 
BC 1.00 1.00 1.00 1.00 


Another factor for the reduced mobility of Sb is the high content of As in the soil, 
which plays a limiting role at 60 - 230 ug / g As (Kabata-Pendias, 2011), and in the 
considered profile the values for arsenic are in the range 137 - 334 ug / g. The mechanism of 
upward movement of arsenic, antimony, selenium and tellurium, which in acidic environment 
are strongly adsorbed by soil colloids as negatively charged oxygen radicals (Sadiq, 1995, 
Denys, 2009, Missen, 2020, Cuvardic, 2003, Plant et al., 2003) is debatable. The measured 
acidity of pH 3.8 creates conditions for the destruction of soil colloids, which in the season of 
wetting, attracting with them the adsorbed ions rise and accumulate in the surface horizon. 
Confirmation of this thesis is the accumulation of physical clay in this horizon, which is more 
than the one below (Table 2), as well as the presence of pyrite and manganese stains 
(Stoykova and Teoharov, 2021). 

The gravelly-sandy composition of the soil-forming sediments and the primitive 
structure of the polluted Regosols presented in profile 8 are a prerequisite for facilitating the 
rise of seasonal groundwater at the foot of the mountain slope and thus covering the entire 
soil profile from glaciation processes. On the same soil-forming materials in the near profile 
7 a cambic soil-forming process takes place, which leads to a concentration of antimony only 
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in the cambic partial carbon horizon (CC = 1.1), again as a manifestation of its metallic 
properties, in association with iron, which is also maximum there (Table 2, Table 3, Figure 
1). 


PROFILE 3 PROFILE 6 


MAsp2k @B1 W1C S2Cg MAsp2k EiaCf CIl SCiilg 


1.50 2.00 2.50 3.00 
PROFILE 7 PROFILE 8 
Ak &Bw(g) GC MACIg MACII| BACIVg 
SE SE 
Sho Gxreebarietarieeereeaers 
Patatatatatetatatatl 
a 
rt 
MIN peeneienunoneueneeennaseuaiy 
0.00 0.50 1.00 1.50 2.00 2.50 3.00 0.00 0.50 1.00 1.50 2.00 2.50 3.00 
PROFILE 11 PROFILE 13 C 
Asp3k ACI CII BA &Biw MBC 


Figure 1. Bar charts of the Concentration coefficients (CC) of Te, Sb, Se and As in 
this study. 

A general trend is shown by the four elements in their redistribution in Rhodic 
Cambisol, presented in the uncontaminated control profile 13c, as the concentration 
coefficients (CC) indicate their downward movement along the profile and their secondary 
accumulation in the BCk horizon. Their downward movement is facilitated by the alkaline 
soil reaction (avg. pH 7.8), in which negatively charged oxygen radicals of metalloids are 
more mobile (Sadiq, 1995, Denys, 2009, Missen, 2020, Cuvardic, 2003, Plant et al. 2003). In 
the surface A-horizon their concentration is the lowest, with coefficients CCAs = 0.49, CCSe 
= 0.49, CFSb = 0.66, CCTe = 0.34. Antimony is again an exception, as it relatively 
accumulates in the A horizon, relative to the intermediate Byw horizon (CC = 0.58). This 
pronounced general tendency to impoverish the upper soil horizons in the unpolluted profile 
is in contrast to the erratic distribution of the four elements in the technogenic contaminated 
profiles (Figure 1). Technogenesis has influenced in such a way that in all contaminated 
profiles no microorganisms were identified, which are known to be an important factor in the 
processes of conversion of metalloids in the soil (Fatoki, 1997, Kabata-Pendias, 2011). 
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Microorganisms are known - fungi and bacteria, that biologically transform metalloids 
through reactions such as oxidation, reduction and methylation, leading to the formation of 
soluble or volatile organic compounds (Fatoki, 1997, Plant et al., 2003, Bentley and 
Chasteen, 2002). Methylation is also known to be a common mechanism by which cells 
metabolize metalloids and is key to the synthesis of volatile compounds leading to depletion 
of soil elements (Fatoki, 1997, Doran and Alexander, 1977, Missen, 2002, Bentley and 
Chasteen, 2002, Denys, 2009). In this sense, it could be reasonably argued that soil 
microorganisms, in the control example, but in the whole study, through their absence or 
presence are a leading factor in the distribution of arsenic, selenium, antimony and tellurium 
in soil profiles. Possible future remediation activities could be based on the pronounced 
biological sensitivity of these elements. 

In profile 3 (Technosols, cambic), contaminated with agglomeration dust, rich in 
heavy metals in the form of technogenic salts (sulfates and carbonates) and Fe-Mn oxides and 
hydroxides and pH 8.4, only antimony has a CC above 1, i.e. it accumulates relatively in the 
surface horizon, and this is its highest concentration coefficient (CC) for all profiles — 2.12. 
Under these conditions, antimony is a companion of Mn, Fe, Pb, Mo, Cu and others 
technogenic pollutants (Stoykova, 2021), attached to Fe-Mn oxide and hydroxide structures, 
exhibiting their pronounced metallic properties, compared to other elements (Stefanova, 
2005, Denys, 2009). 

Table 4. Correlation coefficients (r) between test elements, reference macroelements 
and soil fractions. 


coarse 
Re Wn Se Re te ee dis 
medium sand 
sand 
Fe I 0.48 0.12 0.29 -0.40 -0.36 -0.42 
Mn -0.28 -0.49 
S 1.00 0.20 -0.18 -0.30 -0.48 
As 1.00 -0.21 -0.18 
Se 0.02 0.06 -0.09 
Sb 0.14 -0.15 -0.46 
Te 1.00 0.24 0.06 -0.11 -0.21 


All the presented features are reflected in the derived correlation coefficients between 
the discussed p-elements and some benchmark macroelements - iron, manganese and sulfur. 
Confirmation of the strong attachment of antimony to iron and manganese are the calculated 
high values of the correlation coefficient - r (Sb-Fe) = 0.85, r (Sb-Mn) = 0.62 (Table 4). Man- 
made pollution with agglomeration dust also has a share in the formed large positive 
dependence with iron. Statistically significant, but weaker, is the relationship between arsenic 
and manganese (r (As-Mn) = 0.56). This dependence, as well as the Sb-Mn correlation, are a 
reflection of the influence of coarse-grained soil-forming materials on poluted Regosols 
containing ore. The same is reflected in the statistically significant positive indicators for the 
relationship between the three elements (Mn, As, Sb) and the large soil fraction. The highest 
calculated correlation coefficient is between selenium and tellurium - r (Se-Te) = 0.92. This 
indicator, together with statistically significant values for a positive correlation between the 
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two elements and sulfur, illustrates the primary influence of the halogen geochemical nature 
of these elements, as well as the influence of proluvial soil-forming materials rich in syenite 
and ongoing overwetting processes. Arsenic and antimony are also mobilized in these 
processes, but, although positive, their indicators for correlation with sulfur are not 
statistically significant. 

Conclusion 

After interpreting analytical data from 20 soil samples from Regosols and Cambisols, 
modified in Technosols, extremely high values for arsenic, antimony, selenium and tellurium 
content were determined compared to selected reference values - Te (56x) >> As (25x)> Se 
(24x) >> Sb (10x). The determined average contents are (in ug / g) As — 170,66 (409.34 + 
24.63); Se — 10.56 (14.25 + 6.15); Sb — 6.66 (13.93 + 3.13) and Te — 4.49 (7.03 + 2.14). It 
was specified that the increased contents in these soils are the result of inheritance from 
proluvial soil-forming rocks, fed by two provinces - with predominant syenite or mineral-rich 
composition, and soil-forming and technogenic processes are subordinate. The natural soil- 
forming factors that affect the distribution of these related elements are the overwetting and 
the state of the soil microbiome. Contamination with agglomeration dust is the cause of 
antimony accumulation in the spolic horizon. Antimony, the only one, shows a deviation 
from the common geochemical trends for the four elements, as reflected in the statistically 
significant large positive dependence with iron - r (Sb-Fe) = 0,85. The determined very high 
natural content of the metalloid arsenic in the studied soils, including the arable ones, 
compared to the normatively determined maximum values in Bulgaria, is of public 
importance. The content of arsenic in Technosols is on average about 4 times above the 
maximum permissible concentration of metalloid in industrial soils (40 ug / g), and the 
arsenic in the control profile exceeds 1.5 times the established norm for the maximum 
permissible concentration in arable land (25 wg / g ). Due to the high biological sensitivity of 
the four elements, possible future remediation activities could be based on microbiological 
and organic treatment in situ, following the necessary preliminary laboratory experiments. 
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Abstract 

Due to the detrimental effect of water erosion, erosion data must be gathered for 
sustenance of soil structure and utilization. The aim of the study was to determine the 
erodibility status of the some soils in Akoko-Edo Local Government Area of Edo state, 
Nigeria. Soil samples were taken from representative pedons dug from the three communities 
namely Ikpeshi, Unem-Nekhua and Ososo. The samples were subjected to routine analysis 
and results were used to calculate the erodibility indices of Clay Dispersion Ration (CDR), 
Clay Disperion Index (CDI), Clay Flocculation Index (CFI) and Bouyoucous Erodibility 
Index (ElIrpom). The data derived were subjected to statistical analysis, to determine mean 
separation at 5% probability, coefficient of variation and standard error. The results showed 
that the soils were dominated by the sand fraction, mainly sandy loam with clay increment 
down the horizons becoming loamy sand and sandy clay loam. Soil pH ranged from an 
average of 5.36 — 5.62 across the pedons. Electrical Conductivity (EC) mean across the 
pedons ranged from 23.88 wS/cm (0.02388 dS/m) - 36.11 wS/cm (0.0361 dS/m) which posed 
no salinity threat. Organic matter was low with average of 9.28 — 14.85 g/kg across the 
pedons. Across the horizons, the CDR ranged from 53.7 - 62.3%, CDI (42.5 - 59.2%), CFI 
(34.52 - 43.1%) and Elrom (1.92 - 3.93). Clay had highly significant negative correlation (r= 
-0.7235) with CDI and Elrgom (r = -0.9730). Sand had positive correlation with CDI (r = 
0.7597) and Elgom (r = 0.9130). CDR also had a positive correlation (r = 0.8306) with CDI. 
CDI had a positive correlation (r = 0.6945) with Elrom. Variation for all index was low, 
ranging from 0.1-3. The result shows that the soils are erodible and sustainable and 
conservational agriculture practice should be carried out. 


Keywords: Akoko-Edo, erodibility, Edo state, Nigeria, status 


Introduction 

The need to gather information on soil conditions cannot be overemphasized, as there 
is an up rise in the use of land indiscriminately and also misappropriate land use systems. Soil 
erosion is dependent on erosivity by rainfall and the erodibility which is the resistance of the 
soil to destruction by external forces. The erodibility of any soils depends on the soil type, 
texture, level of aggregation, infiltration rate etc. Areas like Ogbogoro community in Bayelsa 
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State, devoid of soil erodibility data are currently in distress, as the area is ravaged by severe 
erosion. 

As stated by Nyakatawa et al. (2001), soil erosion is a major environmental problem 
worldwide. 

More than 56 % of land degradation is caused by soil erosion, raising a global 
concern on land productivity (Elirehema, 2001). Water erosion conveys soils, its nutrients, 
pesticides and other harmful chemicals into nearby rivers, streams and eventually polluting 
groundwater. Eroded soils are deposited in water bodies leading to pollution and siltation, 
which causes reduction of water quantity, quality, siltation and drying up of water bodies. 
Thus, aquatic life is threatened and eventually eliminated, while land productivity is 
negatively affected. Soil erosion is one of the most serious forms of land degradation in the 
world (Sohan and Lal, 2001). 

Another problem associated with erosion include economic, political, social and 
environmental downturn which result from direct or indirect damages to the soil; which does 
not only cause severe land degradation and low soil productivity but also threatens the 
general health of society, sustainable development and habitation of rural areas, with the 
dwellers in particular (Tang, 2004; Zheng et al., 2004; Jing et al., 2005). 

Consequently, erodibility data must be provided, in order to enable agriculturists, 
engineers and other land users to identify areas prone to erosion. This will enable them 
practice sustainable agriculture, and suitable land use systems, which will help control the 
soil erosion and its hazards. Hence, the study was undertaken to assess the current erodibility 
status of some soils in Akoko-Edo Local Government Area of Edo State, Nigeria. 


Materials and Methods 

Study Area 

The study sites are situated in Akoko-Edo local government area of Edo state, which 
has Ogori /Mangogo, Okehi, Adavi and Okene local government areas all in Kogi state to the 
North; Akoko South East, Akoko North East and Ose local government areas in Ondo state to 
the West; to the South: Owan East and Etsako west; and to the East with Etsako East.. It 
occupies an area of 1,371 km? (137,100 hectares) and lies between longitude 6° 06' 0" E and 
latitude 7° 17'0" N. 

The three areas (Ikpeshi - N 07.18674°, E 006.15456°, Unem Nekhua - N 07.33174° E 
006.08360° and Ososo - N 07.40648°, E 006.25318°) were considered. 

Climate 

Akoko Edo Local Government area is characterized by tropical climate with an annual 
average rainfall amount of 1200 - 1500 mm, mean annual temperatures range of 27° C to 32° C 
and mean annual relative humidity ranging from 30.5 % to 94.0 %. The study area is 
characterized by two distinct seasons: the wet and the dry. The rainfall is at its climax in July and 
August with a short break in mid-August. The dry season begins early November and ends by 
March (Olowojoba et al., 2016). 

Geology 

The area _ lies within the basement complex formation, which consists of various 
minerals like shale, coarse grained granite and granite gneiss etc with some outcrops visible 
around the entire local government area. Physiographically, the project area can be described as 
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situated on gentle plains of low relief in some areas of the local government while other areas 
are quite steep due to the presence of high hills visible in the area. (Olowojoba et al., 2016) 


* Study 


areas 


Fig. 1. Map of Akoko Edo Loca Government Area showing study areas 


Land use 

The major land use system in the study area comprises of livestock farming, arable 
farming, and firewood/ timber exploitation. Arable farming consists of subsistence farming which 
is characterized by intensive and continuous cultivation. The main crops grown in the area are 
cassava, maize, yam, cacao, oil palm, cashew and some other economic trees for timber 
production. 

Vegetation 

The study area is in a derived savannah zone characterized with scattered trees such as 
Irvingia gabonensis, Parkia biglobosa (Jacq) R.Br (African Locust beans), Citrus Sinensis L 
(Orange) and grasses such as Cymbopogon citratus L. (Lemon grass), Pennisetum purpureum, 
Sida acuta Burm F. etc, with some visible rock outcrops occupying some large parts of the local 
government area. 

Soil sampling 

Three profile pits were dug, representing three mapping units/area including: Ikpeshi, 
Unem-Nekhua and Ososo. Soil samples were taken from each distinctive horizons of the soil 
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profiles according to Schoeneberger et al. (2012) guidelines. A total of seventeen (17) soil 
samples were acquired from the three pedons, and were subjected to routine and special 
laboratory analysis. 

Laboratory analysis 

Particle size distribution was determined by Bouyoucos hydrometer method (Gee and 
Or, 2002). Soil pH and Electrical conductivity was determined electrometrically using the 
glass electrode pH and EC meter in a solid-liquid (water) ratio of 1:2.5 (Hendershot et al., 
1993). Organic carbons were determined using Walkley-Black wet oxidation methods 
(Walkley and Black, 1934). The organic matters were computed by multiplying the value of 
the organic carbon by a value of 2.0 by Douglas W. Pribyl (2010). Exchangeable bases were 
determined by the neutral ammonium acetate procedure buffered at pH 7.0 (Thomas, 1982). 
Exchangeable acidity was determined by a method described by McLean (1982). 
Soil bulk density was determined by the core method (Grossman and Reinsch, 2002). Total 
porosity (Po) was obtained from bulk density (fp) values with an assumed particle density 
value of (€s) 2.65 gcm-3 as follows, Porosity (Po) = 100 — (€p/fs) <100/1. 

Erodibility index 

The new and advanced Bouyoucos (1962) erodibility index will be used to analyze the 
erodibility status 


% SAND + % SILT 


Elko => (% CLAY) 


Table 1. ‘ROM’ Scale with regards to soil erodibility category (Roslan & Mazidah, 2002). 


ROM SCALE SOIL ERODIBILITY CATEGORY 
<1.5 Low 
1.5 ~ 4.0 Moderate 
4.0 ~ 8.0 High 
8.0 ~ 12.0 Very high 
> 12.0 Critical 


The clay-dispersion indices were calculated as follows; 

Clay dispersion ratio (CDR) = {(% silt + clay (H20)}/ {% silt + clay (Calgon)} x 100 

Clay dispersion index (CDI) = {[% clay (H20)] / [% clay (Calgon)]} < 100 

Clay flocculation index (CFI) = [% clay (Calgon)] - [% clay (H20)]/ [% clay (Calgon)] < 100 
The higher the CDR and CDI the more the ability of the soil to disperse while the higher the 
CFI, the better aggregated the soil. Clay dispersion ratio was used to determine the erodibility 
of the soils in which greater than 15% are erodible and less than 15 % are not erodible 
(Middleton, 1930) 


Statistical Analyses 

The data generated were evaluated using the Coefficient of variations as designated by 
Wilding et al. (1994) to determine the gradation of disparity among horizons of each pedon. 
GenStat 12" Edition application were summarized using mean, mean separation (Duncan 
Multiple Range Test) at 5% level of probability. The correlation matrix was used to 
determine the association between some selected properties of the soil. 
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Results and Discussions 
Soil Particle Distribution 


Table 2. Some soil physical properties of the study area 


. pilin Clay Silt Sand Textural Class ue poy 
Horizons (cm) Density 
+. g/kg ———_ >» g/cm? %o 
Ikpeshi Soils 
N 07.18674°, E 006.15456° 
Ap 0-24 92.70° 100.00° 807.30° Loamy sand 1.457° 44.914 
A 24-41 104.30° 75.0° 820.70°  Loamy sand 1.743° 34.34° 
Btl 41-76 250° 58.00° 717.00° Sandy clay loam 1.633° 38.49° 
Bt2 76-111 248.33* 60.67° 691.00° Sandy clay loam 1.787° 32.45° 
Bt3 111-147 238.0° 60.00° 702.00° — Sandy clay loam 1.743° 34.34° 
Mean 181.66 70.73 747.6 1.6726 36.91 


Unem-Nekhua 


N 07.33174° E 006.08360 
A 0-14 98.00° 65.00° 837.00° Loamy sand 1.533° 42.26° 
ABI 14-40 105.00° 46.00° 849.00° —_Loamy sand 1.453° 45.28' 
AB2 40 - 64 108.70° 69.30° 822.00° Loamy sand 1.563° 41.134 
BA 64 - 82 90.00' 41.00' 869.00° Sand 1,587° 39.67° 
Btlh 82-118 135.00° 56.00° 809.00° Sandy loam 1.673' 36.98° 
Bt2h 188 — 193 189.30° 51.00° 759.70' Sandy loam 1.653° 34.67° 

Mean 121 54.71 824.28 1.577 39.998 


Ososo soils 


N 07.40648°, E 006.25318° 
0-26 

A 120.00° 95.00" ~— 785.908 Sandy oam 1.283° —-51.73° 
Bwl 26 - 49 205.00° 89.00° 706.00° Sandy clay loam 1.567° 39.62° 
Bw2 49 - 83 238.00° 80.00° 682.00° Sandy clay loam 1.683° 36.63" 
Bw3 83-114 288.00° 75.00° 637.00° Sandy clay loam 1.587° 40.33° 
Bw4 114 - 144 220.70° 183.31° 596.00' — Sandy clay loam 1.607° 39.25" 
Bw5 144 - 187 285.00° 70.00' 645.00° — Sandy clay loam 1.423° 46.25° 

Mean 226.11 98.71 675.16 1.525 42.3 


Mean value(s) with the same letters(s) in the column are not significantly different from one another at 5% level 
of probability in each location. 


The result from table 2 shows that sand fraction dominated all the horizons of the 
Ikpeshi soils, followed by clay, and then silt. The soils had a mean of 747.6 g/kg sand, 181.66 
g/kg clay and 70.73 g/kg. The sand content reduced down the horizon while the clay fraction 
increased downwards. The soil texture was loamy sand at the first two horizons and sandy 
clay loam for the remaining three horizons. The sandy and coarse nature of Ikpeshi soils is 
greatly related to the parent material: basement complex formation, which consists of various 
minerals like shale, coarse grained granite and granite gneiss (Weil and Brady, 2016). The 
increase in clay down the horizon indicates the process of illuviation Niu et al. (2015). In 
Ikpeshi soils, the bulk density ranged from 1.457 — 1.787 g/cm? with an average of 1.67 
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g/cm*. The bulk density of the Ikpeshi soil was high indicating compaction, in exception of 
the first horizon which had been ploughed. The bulk density values affected the porosity of 
the soils, the horizon with low bulk density had lower porosity of 44.91% while the highest 
bulk density of 1.787 g/cm3 registered a porosity of 32.45%; lower bulk density begot higher 
porosity while higher bulk density begot lower pore spaces. 


Table 3. Selected Chemical Properties of the Pedons 


Depth 


Horizons (enn) pH EC Org.C Org.M CEC ECEC 
pwS/em « glkkg —y,»Y  <+— cMol/kg __» 
Ikpeshi Soils 
N 07.18674°, E 006.15456° 
Ap 0-24 4.43° 20.53° 11.93* 20.56° 2.36" 3.16° 
A 24-41 5.73° 24.63° 4.434 7.63° 2.80° 3.43° 
Btl 41-76 6.10" 52.80" 10.90° 18.79° 2,444 2.84° 
Bt2 76-111 5.96° 31.73° 7.93° 13.67° 2.61° 3.04° 
BB 111 - 147 5.90° 24.70° 7.90° 13.61° 3.00° 3.60° 
Mean 5.62 30.88 8.6 14.85 2.64 3.21 
Unem-Nekhua 
N 07.33174° E 006.08360 
A 0-14 4.80° 17.20° 2.50° 4.31° 222° 3.03° 
ABI 14-40 5.60° 39.3" 6.30° 10.86° 2.44° 3.24* 
AB2 40 - 64 5.20° 20.30° 6.208 10.68° 232° 2.92° 
BA 64 - 82 4.30' 6.50f 1.70° 2.93! 3.08° 3.22° 
Btlh 82 - 118 6.60" 24.40° 6.40° 11.03° 2.26° 2.66' 
Bt2h 188 — 193 5.66° 35.63” 9.23" 15.91" 2.33° 2.92° 
Mean 5.36 23.88 5.38 9.28 2.44 3 
Ososo soils 
N 07.40648°, E 006.25318° 
A 0-26 5.50° 49.50° 7.90° 13.61° 2.50° 3.30! 
Bwl 26 - 49 5.50° 40.10b 8.00° 13.79" 3.16° 3.96" 
Bw2 49 - 83 5.10° 29.30° 7.808 13.449 278" 3.58° 
Bw3 83 - 114 5.60° 37.40° 6.90° 11.89° 3.34" 3.94° 
Bw4 114-144 5.23° 32.00° 7.86° 13.55° 3.08° 3.84° 
Bw5 144 - 187 5.30° 28.40' 6.40' 11.03° 3.08¢ 3.68° 
Mean 5.37 36.11 TAT 12.88 2.99 3.71 


EC: Electrical Conductivity, OC-Organic Carbon, OM — Organic Matter, CEC- Cation Exchange 
Capacity, ECEC — Effective Cation Exchange Capacity 


In Unem-Nekhua soils, the sand fraction also dominated the horizons with an average 
of 824.28 g/kg, followed by clay with a mean of 121 g/kg and silt 54.71 g/kg. The texture 
ranged from loamy sand to sand to sandy loam indicating a reduction of the sand fraction and 
clay increment down the horizon. I[lluviation also took place in the soils. The bulk density 
increased down the profile with an average of 1.577 g/em* and an average porosity of 39.99% 

Ososo soils recorded the highest clay content with an average of 226.11 g/kg and 
lowest sand fraction of 675.16 g/kg, the soils were majorly sandy clay loam apart from the 
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first horizon. The soils had the lowest average bulk density of 1.53 g/cm’ and highest porosity 
of 42.3%. 

Some chemical properties of the study areas 

Ikpeshi soils had a pH of 4.43 — 6.10 which is strongly acidic to slightly acidic. The 
high acidity could be linked to the leaching of basic cations caused by the high rainfall in the 
region Niu et al., (2015). The electrical conductivity ranged from 20.53 — 52.80 uS/cm with a 
mean of 30.88 uS/cm (0.03088 dS/m). This value according to the Ganjegunte et al. (2018) 
places the soil less than four (<4), which is non saline, indicating that there is no obstruction 
to the soil structure and is poses no threat to seedlings and salt sensitive crops. Organic 
carbon and organic matter low to moderate with a range of 4.43 — 11.93 g/kg and 7.63 — 
20.56 g/kg. The low effective cation exchange capacity ranged from 2.84 — 3.60. This 
concurred with the findings of Gailyson and David (2013). The low CEC and ECEC 
indicated that the pedons have low activity clays and are highly weathered according to the 
finding of Soil Survey Staff (2014). The highly weathered soils are more susceptible to 
erosion due to the formation of coarse particle which encourages poor aggregation. 

Unem-Nekhua soils had pH with a mean of 5.36, EC of 23.88 wS/cm (0.02388 dS/m), 
Organic Carbon and Organic matter of 5.38 and 9.28 g/kg, Cation Exchange Capacity and 
Effective Cation Exchange Capacity of 2.44 and 3 cMol/kg while Ososo soils recorded a 
mean pH of 5.37, EC of 36.11 wS/cm (0.0361 dS/m), OC and OM of 7.47 and 12.88 g/kg; 
with CEC and ECEC of 2.99 and 3.71 cMol/kg. The interpretations for Unem-Nekhua and 
Ososo were all similar to that of the Ikpeshi soils. 

From the result in table 4, Ikpeshi soils had a Clay Dispersion Ratio (CDR) ranging from 
44 — 67% with a mean of 53.5%. Variation and standard error amongst the generic horizons were 
low with 0.1 coefficient of variation (CV) with standard error (SE) of 0.0577. The Clay 
Dispersion Index (CDI) ranged from 40.4 — 52% with an average of 44.2; the coefficient of 
variation and standard error were 0.1 and 0.0577. The Clay Flocculation Index (CFI) ranged from 
30.7 — 39.9% with a mean of 34.52, Coefficient of Variation of 0.2 and Standard error of 0.0577. 
The Bouyoucos Erodibility Index (Elgom) ranged from 1.52 — 4.88 with a mean of 2.80, CV of 
0.9 and SE of 0.02633. The CDR, CDI, CFI and Elgom registered low variations amongst the 
different horizons. 

The CDR of the Ikpeshi soils was far above 15% in all the horizons indicating that they are 
weak soil structure and are erodible. The CFI which is another factor used to determine the 
resistance of soils to dispersion, was also low indicating low resistance of the soils to external 
erodible forces. The CFI connotes the ability of the soils to aggregate; the Ikpeshi soils recorded 
low CFI, thus weak aggregation. The Elrom on the other hand ranged from high to moderate 
down the horizons. 

The CDR in Unem-Nekhua soils ranged from 50.5 — 77% with a mean of 62.3%, CDI (50 
— 66.2%, mean — 59.2), CFI (33.5 — 40.89% with mean of 36.4) and Elrom had a range of 2.3 — 
5.1 with a mean of 3.90. The coefficient of variation and standard error were 0.1 for CDR, CDI, 
CFI respectively, 1.5 for Elgom and 0.0473 for CDR, 0.0577 for CDI, 0.0410 for CFI, & 0.0577 
for Elrom. The table 4 shows that CDR and CDI were above 15% indicating that the soils are 
erodible, also the CFI for the different horizons were low, showing weak aggregation. The Elrom 
indicated that the erodibility was high to moderate. 
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Table 4. Erodibility indices of the study areas 


ee CDR CDI CFI igen 
Horizons 
(cm) ——_—_————_- % — > 
Ikpeshi Soils 
N 07.18674°, E 006.15456° 
Ap 0-24 67° ele 39.9° 4.88° 
A 24-41 44° 40.4a 30.7a 431d 
Btl 41-76 s7a° 458° 35.7¢ 1.72° 
Bt2 76-111 S14" 42.5° 33.3° 152° 
Bt3 111-147 49° 40.3° 33° 1.6 
Mean 53.7 44.2 34.52 2.80 
CV 0.1 0.1 0.2 0.9 
SE 0.0577 0.0577 0.0577 0.02633 
Unem-Nekhua 
N 07.33174° E 006.08360 
A 0-14 ae 66.2' 34.01° 4.6° 
ABI 14-40 50.52° 50° 40.89' 4:38 
AB2 40 - 64 56.8° 55.1° 39.02° 4.1° 
BA 64 - 82 67.9° 65.8° 33.5" a 
Btlh 82-118 61.41% 60.7! 34.8° 32° 
Bt2h 188 — 193 60° 576° 36° 2.3° 
Mean 62.3 59.2 36.4 3.93 
CV 0.1 0.1 0.1 is 
SE 0.0473 0.0577 0.0410 0.0577 
Ososo soils 
N 07.40648°, E 006.25318° 
A 0 — 26 57° 50.5° 45.05° aT" 
Bwl 26 - 49 56.2° 46.9° 42.69 1.9° 
Bw2 49 - 83 56° 423° Ay 1.6° 
Bw3 83-114 56.01° 40.8° 41.7° 1.28 
Bw4 114-144 50.3° 40.8° 41.6° 1.8° 
Bw5 144 - 187 493° 33.7" 45.8! is 
Mean 54.1 42.5 43.1 1.92 
CV 0.1 0.1 0.1 3.0 
SE 0.0528 0.0577 0.0528 0.0577 


CV — Coefficient of Variation 
Elrom - Bouyoucos Erodibility Index 


SE — Standard Error 


CDR — Clay Dispersion Ration 
CDI — Clay Dispersion Index 


CFI — Clay Flocculation Index 


In Ososo soils, the CDR ranged from 49.3 — 57% with a mean of 54.1% , CV (0.1), and 
SE of 0.0528. The CDI was 33.7-50.5% with a mean of 42.5%, CV of 0.1 and SE of 0.0577. 
CFI ranged from 41.6 - 45.8% with a mean of 43.1%, CV of 0.1 and SE of 0.0528. Elrom 
ranged 1.3-3.7 with a mean value of 1.92, CV of 3.0 and SE of 0.0577. The data shows that, 
like the other soils, the Ososo soils are also erodible. Igwe and Udegbunam (2008), also 
reported that DR, CDI, and CFI were high when compared to the soils of Southern Nigeria. 
The result when correlated to the findings of Oguike and Mbagwu (2009) in soils of Southern 
Nigeria shows that DR and CDI were high while CFI was low. 
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Table 5 shows the correlation matric between some erodibility indices and some 
selected soil physical properties. Clay had highly significant negative correlation (r= -0.7235) 
with CDI and Elpom (1 = -0.9730). Sand had positive correlation with CDI (r = 0.7597) and 
Elrom (rt = 0.9130). CDR also had a positive correlation (r = 0.8306) with CDI. CDI had a 
positive correlation (r = 0.6945) with Elrom. 


Table 5. Correlation matric of erodibility indices and some selected soil physical properties 


CLAY SILT SAND CDR CDI CFI Elrom 
——___ gkkg ————> +———_ % ———> 
CLAY 1 
SILT 0.1507 1 
SAND -0.9227** -0.5199 1 
CDR -0.4585 -0.1878 0.4664 1 
CDI -0.7235** -0.3490 0.7597** 0.8306** 1 
CFI 0.2917 0.4248 -0.4167 -0.1374 -0.3473 1 
EIrgom -0.9730** -0.1858 0.9130** 0.4918 0.6945** -0.2713 1 
Conclusion 


The results of the study areas show that the soils were dominated by the sand fraction 
which could be attributed to the parent material, the pH was acidic to slightly acidic; the soil 
of the three areas recorded low electrical conductivity indicated no threat of salinity and 
obstruction. The organic matter and cation exchange capacity was also low. The erodibility 
indices used stated that the soils are erodible. This could be related to the low organic matter 
content and lower clay fraction, which could serve as soil binding agents and improve soil 
structure and aggregation. The destructive power of erosion on agriculture can never be 
overemphasized, hence conservation methods such as cover cropping, contour cropping, 
slope terracing, and conservational tillage methods etc should be practiced to control and 
limit soil erosion. 
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Abstract 

The soils from the Dolna Banya valley of Sofia district are formed under forest 
meadow, bush and tree vegetation. It has influenced their texture, natural fertility and the 
reaction of the soil solution. As they are located in a semi-mountainous terrain, surface 
horizons are subject to erosion processes and podzolization. 

The aim of the research is to establish the changes in the profile of Rhodic Luvisols 
which occurred as a result of natural influences and the type of production activity on them. 

The research has found that for a long time there have been changes in the surface soil 
horizons of all studied differences of Rhodic Luvisols. The process of podzolization 
continues to develop and there is an increase in the content of silt in the horizons below the 
humus horizon. In the case of Fluvisols, an increase in the fraction of coarse and fine sand is 
reported. In the case of moderately eroded Planosols land use has an impact on the soil 
texture. In the deeper horizons, when planted with raspberries, the content of silt and physical 
clay increases compared to the same soil difference in the area with the application of three- 
field crop rotation, compacted with cereals using a shallow mouldboard ploughing . 


Keywords: Fluvisols, Luvisols, Planosols, Sofia district, process of 
pseudopodzolization, soil texture, three-field crop rotation 


Introduction 

From the fact that the properties of the soils and the degree of their fertility depend a 
lot on the type and age of the parent rocks, on the relief, exposure, latitude and altitude, as 
well as on the climatic conditions that activate or suppress more or less biological activity, 
the state and biological potential of soil differences are determined (Teoharov 
2004; Blum, 2007; Kinoshita et al., 2017). This feature of the soil formation process is an 
objective basis for the action of the third law of agriculture: "for the unity and 
interdependence of growth and development of plant associations and their habitat” from 
which arises the objective need for an ecological approach in organizing and managing 
production in crop production (Krastanov, 2006; Ugarte et al., 2016). 

The soils from the Dolna Banya valley are formed under forest meadow, bush and tree 
vegetation, which has affected their texture, natural fertility and the reaction of the soil 
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solution (Koynov et al., 1972; Penkov et al., 1985; Ninov et al., 1985; Boneva et al., 2005; 
IUSS Working Group, 2015; Teoharov et al., 2019). 

As they are located on a semi-mountainous terrain characterized by rough terrain with 
ravines and gullies that have slope 8-10° and the soil surface horizons are subject to erosion 
processes and pseudopodzolization. The climate is transition continental region of the 
European continental climate and is temperate continental influenced by the river valley of 
Maritsa river and its tributaries (Boneva et al., 2005; Ninov, 2007). 

Economic activity in the past was determined by the cultivation of cereals, fiber crops, 
potatoes and mainly berry crops, natural and artificial pastures and meadows. 

The aim of the research is to establish the changes in the profiles of Luvisols, 
Fluvisols and Planosols which occurred as a result of natural influences and the type of land 
use on them. 


Materials and Methods 

The subject of the study are the following differences of soils: 

Rhodic Luvisols, slightly eroded. 

They are distributed in the form of individual spots in the southeastern and 
northeastern part of the land. They occupy an area of 3875 da - which represents 2.39% of the 
total area. They are formed on non-carbonate weathering products under the influence of 
forest vegetation. As a result of the recent slight erosion, part of the light surface horizon has 
been removed, due to which the profile is somewhat shortened. The profile is clearly 
differentiated. Most of the area of the considered soil differences is included in the arable 
fund of the farm units. An idea of the morphological structure of these soils is given by the 
description of profiles made one kilometer northeast of the village of Pchelin (Yolevski, M. et 
al. 1974). 

AB fallow land 0-20 yellowish-brown (7.5YR/5/4 in the dry state), fresh, 
compacted, clayey-sandy, granular and crumbly structure, does not carbonate from 
hydrochloric acid, transition noticeable. 

B,;, 20-36 reddish-brown (SYR/3.5/ 6 in dry state), rusty-reddish spots, slightly 
moistened, dense, medium sandy-clayey, prismatic-lumpy structure, not carbonated by 
hydrochloric acid, transition noticeable. 

Bo; 36- 54 reddish-brown (5YR/3/6 in the dry state), rusty-reddish spots, moist, dense, 
medium sandy-clayey, prismatic-lumpy structure, not carbonated from hydrochloric acid, 
transition noticeable. 

BC 54-70 variegated (SYR/4/4 and 5YR/6/4 in the dry state), moist, strongly 
compacted, medium sandy-clay, lighter than the above, unstructured, does not carbonate from 
hydrochloric acid. 

The humus horizon has a thickness of 18-20 cm. 

The soil texture in the surface horizon is clayey-sandy. The amount of physical clay is 
17.2%. The soil texture in the illuvial horizon is medium sandy-clayey 39.7 to 42.1%. The 
amount of silt in the surface horizon is low, while in the illuvial horizon it is several times 
higher - 32.7 - 36.1%. The predominant fractions are those of silt, medium and fine sand. 
This soil texture determines poor water permeability, good water retention and unfavorable 
air regime. According to the amount of organic matter in the surface horizon, the soil is with 
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low humus content -0.83%.The have a low reserve of total nitrogen and total 
phosphorus. Carbonates are not contained in the depth of the entire soil profile. The soil 
reaction has a strongly acidic pH in H2O of 5 to 5.3. These soils have a number of 
unfavorable properties that make them suitable for a limited number of crops - very low 
humus content, sharp textural differentiation, strongly acidic soil reaction (it can be corrected 
by liming). They are suitable for growing some less demanding crops. 

Eutric Planosols, moderately eroded 

They occupy an area of 7688 da, 5.15 % of the total area southeast of the village of 
Dolna Banya.They are distributed mainly in the central flattened part of the territory of the 
farm units. They are formed on Pliocene and Tertiary deposits under the influence of forest 
vegetation. No significant erosion processes are observed in them. Their morphological 


structure can be observed by the description of a soil profile on a plateau in fruit orchard 
located 2 km north of Dolna Banya village (Yolevski, M. et al. 1974). 

AiE 0-29 Light gray-brown, fresh, loose, clayey-sandy, small crumble structure, 
powdery, not carbonated from hydrochloric acid, transition clear. 

AoB(t)g 29-42 light brown, fresh, compact, slightly sandy-clayey, crumble structure, 
partially dusted by pseudopodzolization, dusted with amorphous silica, not carbonated from 
hydrochloric acid, simple transition. 

B,(t)cg 42-75 light brown-reddish, moist, slightly dense, slightly sandy-clayey, 
prismatic, signs of gleying, contains manganese nodules of grains, does not carbonate from 
hydrochloric acid, transition gradual. 

B2(t)eg 75-115 brownish-reddish, moist, slightly dense, heavy sandy-clayey, 
prismatic, contains manganese nodules, does not foam from hydrochloric acid, transition 
noticeable. 

C 115-140 light brownish-white, moist, compact, slightly sandy-clayey, unstructured, 
gravelly, slightly carbonated from hydrochloric acid. 

Morphologically, these soils are characterized by weak (20-25 cm) humus horizon 
with light brown color and dusty from pseudopodzolization and processing small crumbly 
structure. Dusting with amorphous silica is observed here. Beneath it lies a deep about 70-100 
cm alluvial horizon with a reddish dense color, clayey with a prismatic structure and the 
presence of manganese nodules on water-bearing parts, this horizon is less deep and contains 
unweathered rocks. Horizon C is almost intact Pliocene Tertiary deposits with light brown to 
off-white color. The terrain occupied by these soils is flat (with a slight slope) convenient for 
farm machines. 

The soil texture in the surface horizon is clayey-sandy and the amount of physical clay 
is 18.4%. In the alluvial horizon the soil texture is sharply aggravated - the amount of 
physical clay is 51.1-62.8%. 

The amount of silt in the upper two horizons is low while in the range of the illuvial 
horizon it is 45.6-57.4%. The predominant fraction of soil profile depth is that of silt followed 
by the fraction of medium and fine sand. The texture coefficient is 7.8. 

This soil texture causes low water permeability, high water holding capacity and 
unfavorable air regime. The content of humus 0.94% in the surface horizon defines these 
soils as having low humus. In the lower divisions of the soil profile the amount of organic 
matter is insignificant - 0.36-0.5%. Carbonates are not contained in the depth of the entire soil 
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profile. The soil reaction is medium acidic in the surface horizon - pH 5.8 (H2O). With good 
farming techniques some less demanding crops can be grown on these soils. 

Eutric Fluvisols 

They occupy an area of 2919 da - 1.96% of the total area. 

They are formed on Alluvial-Deluvial deposits under the influence of meadow vegetation. 
They have a layered structure characteristic of all alluvial soils. The materials creating the 
individual layers have a light soil texture. Humus is better expressed in the upper part of the 
soil profile. The morphological structure of these soils is given by profile 2: 

A; 0-11 light gray-brown (7.5YR/5/6 in the dry state), fresh, loose, slightly sandy- 
clayey, weakly expressed fine-grained structure, does not carbonate from hydrochloric acid, 
transition noticeable. 

Aq layer 11-24 grayish-brown (7.5YR/5/4 in the dry state), slightly moistened, loose, 
sandy (loose sand), unstructured, does not carbonate from hydrochloric acid, transition 
noticeable. 

Amt layer 24-50 grayish-brown lighter than above (7.5YR/5/6 in the dry state), slightly 
moistened, loose, sandy (loose sand), unstructured, does not carbonate from hydrochloric 
acid. 

The humus horizon has a thickness of 20-30 cm, the amount of organic matter in the 
humus horizon is 2.17-3.19%. The depth of the soil profile is 60-90 cm. Soil texture in the 
surface horizon is clay-sandy, the amount of physical clay is 20.1-28.5%. The predominant 
fractions are coarse particles. This soil texure determines very high water permeability, low 
water holding capacity and favorable air regime. 

The soil reaction is moderately acidic to neutral pH in KCI with value of 5-5.7. 

Leptic Cambisols, eroded, stony 
Occupy an area of 15688 da - 10.5% of the total area of the study area.They are formed on 
weathering products on acid rocks under the influence of forest vegetation and active 
erosion.The profile of these soils is very short and the hard rock is closely laid. They have a 
light soil texture; they contain a high percentage of skeletal material and a small percentage 


of organic matter. Due to the very large slope occupied by these soils they are almost 
economically irrelevant since only individual sections are used. 

Our research is based on the analysis of soil samples from the area representing the 
natural meadow (profile 1) , winter cereal - rye (profile 2) with raspberry plantations (profile 
2) and left in the prologue, after some years of use ( profile 3). 

Researched indicators: 

Physical and texutral parameters: 

- soil texture - by pyrophosphate method (by Kaczynski, 1958) ; 

- structure and waterproof aggregates - by sieve analysis dry and wet sieving by the 
method of NI Savinov (Penkov et al. 1997). 

«Agrochemical indicators: 

-mobile forms: nitrogen - by the method of Bremner and Kiney; phosphorus - by the 
method of P. Ivanov; potassium - P. Ivanov (Penkov et al. 1997); 

-humus - by the Turin method (Kononova, 1963); 

-pH - potentiometrically (in water and in potassium chloride), (Penkov et al. 1997); 

-carbonate content - by the Scheibler method (Penkov et al. 1997). 


Bulgarian Journal of Soil Science® 2022 Volume 7. Issue 1. www.bsss.bg 27 


Kirilov et al., 2022 


An agrotechnical assessment was performed for the condition of the production areas 
subject to the study - soil condition in the humus horizon, presence of plant residues, weed 
density, etc. 

Results and Discussion 

From the comparative analysis of the results of the laboratory determination of the 
soil texture, changes in the fractional composition of all studied soil differences are 
established. For a period of over fifty years, depending on the land use, microrelief and the 
influence of meteorological conditions, the changes in the main soil horizons are more or less 
distinctive. 

Table 1. Soil texture profile 1 and 54 Dystric Planosols 

Horizon and 1.00 - 0.25 - 0.05 - 0.01 - 0.005 - 


a - 36-54* 


ole T 


A,-0-20** 


* Note - the data are from the saeeE of ae in 1974, Yolevski Moe et al. ,,Soil characteristics of the land 
of Dolna Banya — Sofia district, archive of IP ’N. Pushkarov’, Sofia. 
** Note - the data are from the analysis of samples in 2019. 

A comparative analysis of the data on the soil texture of Dystryc Planosols shows that 
despite the long-term difference between the two studies, the fluctuations in the fractional 
composition are insignificant (profile 1) . 

The long-term impact and sporadic production activity on the land plot did not 
significantly affect the soil texture of the soil disctinction. In the surface horizon, the silt and 


physical clay content was almost unchanged for the long period between the two 
analyzes. More significant differences between the fractional composition were found in the 
B, horizon. The content of the silty fraction decreases sharply - by 10%, and the content of 
physical clay decreases by 8%. Probably in the C horizon there is a porous rock mass in 
which parts of the fine fractions of the clay have deposited. 

Of note is the fact that as a result of the absence of land use the humus content 
increases - on average by 0.09%. The inclusion of the area in the production activity was for 
a short time and did not create a risk of surface water erosion. 

When comparing the data available from the soil outline and those of made 
agrochemical analysis indicated that the mobile forms of nitrogen in arable layers is reduced 
to 4-5 mg / kg soil, and that of useful phosphorus 2.5 times. There is an increase in the useful 
forms of potassium in the surface horizon, probably due to the decomposition of straw from 
cultivated cereals. The reaction of the soil solution remains moderately acidic in the surface 
horizon, but in the B horizon it becomes slightly acidic. 
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Table 2. Agrochemical properties profile 1 and profile 54 Dystric Planosols. 


mg/kg mg/ 100 ¢ Humus 


Horizon and Base 


depth Saturation 
(V%) 


5.0 


Bit, 20-36* 


Bot, 36-54* 5.1 


Bot, 36-54** | 5.5 
* Note - the data are from the analysis of samples in 1974, Yolevski, M., et al. ,,Soil characteristics of the land 
of Dolna Banya — Sofia district, archive of IP *N. Pushkarov’, Sofia. 
** Note - the data are from the analysis of samples in 2019. 


During the long period between the two studies, profile 2 revealed significant changes 


in horizons A and B, which are a consequence of the impact of climatic factors and its 
inclusion in the land fund. Due to the application of a crop rotation system of agriculture 
through the alternation of cereals oats and rye and trench flax and potatoes, the unification of 
the soil texture at a depth of up to 42 m is established. 

Table 3. Soil texture profile 2 and profile 4 Eutric Planosols. 


Horizon depth >1 1- 0.25- | 0.05- | 0.01- | 0.005- | <0.001 | ¥ <0.01 
(cm) 0.25 | 0.05 0.01 | 0.005 | 0.001 


Ase a [4a0 [14 [a3 | 162 | 136 | 69 | 32 | 131 | 272 | 
ae [| 0 psa pase [9 pas | 20 [ra] | 


Profile 2 


ARID [3G] 00 pwEl ITM [ssT oOo] wl a 

Abe 2-42 | 536 | 00 [371] 56 | 50 | 57 | 65 | 11 | 23 

Be as | - | oo [s2 fia [asa | ios] 53 | 399 | 358 | 
Profile 2 

ARO TT 00 ps4] 1] 34 [59] 68 | Wa] OY 

ABs [ | 00 [aie] 379 | 21 | 59 | 89 | 134 | 282 

Bg | _f oo foi pss [aus] es] 7 | aa [as | 


* Note - the data are from the analysis of samples in 1974, Yolevski, M., et al. ,,Soil characteristics of the land 
of Dolna Banya — Sofia district, archive of IP ’N. Pushkarov’’, Sofia. 
** Note - the data are from the analysis of samples in 2019 from levels with three crop rotation 


*** Note - the data are from the analysis of samples 2019 perennial plantation - raspberries 

Unlike the initial data, where the process of pseudopodzolization is visibly outlined in 
the three horizons, in the analysis in 2019 there is a visible textural boundary between A and 
B horizon, as the content of silt is 4.5 times higher and physical clay over 3 times. In land use 
with perennial plantations - raspberries, a higher content of silt and physical clay is found in 
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the A,B horizon, as compared to the levels with predominant cultivation of cereals and 
cereal-legume mixtures, the physical clay has a 5% increase. A probable reason is the fact 
that during the creation of the raspberry plantation deep plowing was performed, followed by 
plowing in the rows every year while during the three-field crop rotation the depth of 
cultivation is up to 20 cm. 

The chemical properties data for the two surface horizons outline some trends. The 
reaction of the soil solution of weakly acidic to neutral after long period becomes moderately 
acid (5.5-6.2 in H,O). This change is due to natural acidification processes under the 
influence of water from the adjacent forest vegetation and not from production activity. The 
land section has been artificial meadow, which then becomes a natural grass. The soil 
is depleting of bases in the surface horizon, despite the percentage increases in the depth, 
reaching 78.0%. 

Table 4. Chemical properties Profile Ne 2 and Profile Ne 4 Eutric Planosols 


Hori 4 Carbonates Base mg/ 
ronan Saturation kg mg/100 g Humus 
depth 
(em) 2 Total Yo 
N P.O; K,O 


Profile 4 


=[5 
AiBg 29-0 | 
Big 2-75" | 65 |S [—770_[ 96 | 12 [4a | ona | 
Profile 2 

EOD [3s] 4 55 
ABg2 a" 730 
Big 42-75 ca a a 


AE 0-29 *** 
A>B 29 — 42*** 
Btg 42 -75*** 
* Note - the data are from the analysis of samples in 1974, Yolevski, M., et al. ,,Soil characteristics of the land 
of Dolna Banya — Sofia district, archive of IP ”N. Pushkarov’”, Sofia. 
** The data are from the analysis of soil sample from 2019 from a field with three crops rotation 
*** The data are from the analysis of soil sample from 2019 with permanent crop raspberries 


During the long-term land use of the area it is established that during the raspberry 
plantation the assimilable forms of nitrogen increase, as for the horizons to a depth of 


42 cm it is respectively 5.8 mg/kg of soil. In the absence of balanced fertilization, the mobile 
forms of phosphorus and potassium decrease and it is more significant for the amount of 
P,O;- 4.5 times in the area with three crop rotation and about 3 times in the area occupied by 
raspberries. We assume that when creating the perennial plantation, basic fertilization with 
phosphorus and potassium was performed. 

Pseudopodzolic soils differences, in addition to a pronounced textural differentiation, 
are also characterized by a weak humus horizon and an acid reaction of the soil 
solution (Koynov et al., 1972) . The treatment of these soils should be carried out during time 
of appropriate humidity to avoid some adverse physical and mechanical consequences. The 
arable land can be deepened to 25-30 cm with the application of manure and higher doses of 
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mineral fertilizers. When using mineral fertilization priority should be given to fertilizers 
with physiologically alkaline actions. 


The structure of the Fluvisols (Alluvial-Deluvial) soil is layered, with stratums which 
is formed under the influence of the flowing waters of the Maritsa River. Subsequently, the 
river withdrew to its current bed, and a process of deluvial formation began on the terrain. Of 
note is the fact that the percentage of coarse sand is increasing - by 15.7% and 21.8% 


respectively for Ayjlayer. and A ylayer. horizons. 
Table 5. Soil texture profile Ne 3 and profile Ne 46 Eutric Fluvisols 


Horizon Amount | I- | 0.25- | 0.05- | 0.01- | 0.005- | <0.001 | ¥ <0.01 
depth >1 0.25 | 0.05 | 0.01 | 0.005 | 0.001 
(cm) 


profile Ne 46 


P Auaasor | - | 792 [4] 42 [0s | o2 | 02 | 62 | 160 | 
profile Ne 3 

Aor [218 [00 [30 [39 [45 [55] OA TT | BO 

Aura [2.60 [00 [322 [400 | 37 | 55] 64 | 22 | 241 

PAnzaso [ - [00 [ofa] 9s | 66 | 33 | 99 | 28 


* Note - the data are from the analysis of samples in 1974, Yolevski, M., et al. ,,Soil characteristics of the land of Dolna Banya — Sofia 
district, archive of IP *N. Pushkarov’”, Sofia. 
** Note - the data are from the analysis of soil samples in 2019. 


Based on microrelief, it can be assumed that this change in the soil texture is due to 
the surface water erosion with the adjacent land plot of asphalt road the soil was moved to 
declines forming the ravine of the river. A profile made in an adjacent pine grove gives 


grounds for such an assumption. 
Table 6. Chemical properties profile 3 and profile 46 Eutric Fluvisols 
Carbonates Base 
(CaCO3) Saturation 


* Note - the data are from the analysis of samples in 1974, Yolevski, M., et al. ,,Soil 
characteristics of the land of Dolna Banya — Sofia district, archive of IP ”N. Pushkarov’”, Sofia. 
*** Note - the data are from the analysis of soil samples in 2019. 

The chemical properties are characterized with slightly acidic to neutral pH in soil 
research in cooperative farm Dolna Banya which in 2019 showed a decrease to medium 
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acidic soil solution. Carbonates are washed deep into the profile and base saturation slightly 
decreased as a result of the impact of climatic factors. In the event of the water erosionthe 
humus content is reduced of the Eutric Fluvisols, as for long-term decrease was respectively 
0.39% and 0,06% for the first and second layer of the soil profile. 


Conclusion 

For a long period of time there were changes in the surface soil horizons of all studied 
differences of soils. The process of pseudopodzolization continues to develop in Planosols. It 
is observed an increase in the content of silt in the horizons located below the humus 
horizon. In the case of Eutric Fluvisols (Alluvial-Deluvial), an increase in the fraction of 
coarse and fine sand and a decrease in the humus content are reported. Land use has an 
impact on the soil texture in the case of moderately eroded Dystric Planosols. In the deeper 
horizons, when planted with raspberries, the content of silt and physical clay increases 
compared to the same soil difference in the area with the application of three crop rotation, 
compacted with cereals with a subsided surface. The content of mobile forms of nitrogen 
increases in the area with permanent planting. The useful for plants forms of phosphorus and 
potassium decrease more significantly in the area on which three crop rotation cultivated with 
cereals with a mouldboard ploughing for growing is applied. 
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Abstract 

Understanding the influence of agricultural soil management practices on soil 
processes and the environment is vital for sustainable food production. Therefore, this study 
evaluated the impact of soil temperature and moisture on soil organic carbon mineralisation 
under synthetic (NPK) and organic (biochar) amendments in a maize cropping system in 
Torin, Southern Guinea Savanna of Nigeria. The objectives were to determine the C 
mineralisation and C sequestering potential under each amendment and assess how 
precipitation amounts and soil temperature influence carbon mineralisation. Two soil types 
(Anthrodensic Ustorthent and Typic Ustorthent) were used, and four treatments (control, 
NPK, biochar+NPK, and biochar) replicated thrice were laid out in completely randomised 
design on each soil. Soil respiration was measured in-situ by the alkali absorption method. 
Data collected was subjected to analysis of variance procedure for complete randomized 
design to determine significant differences in treatment means, correlation and regression 
analyses were also carried out to test the relationships amongst experimental factors and 
variable responses. Results showed that evaluated nutrient management options did not 
significantly influence C mineralisation. C mineralisation, however, differed between soils, 
with Typic Ustorthent having higher organic carbon (OC), total CO emission (82.35 mg) and 
mean CO, emission rates (1.18 mg m~ d’). Correlation analysis showed differing 
relationships between soil temperature and C mineralisation under the different treatments 
and soils, with significant positive correlation (r=0.44; p= < 0.05) recorded only under sole 
biochar treatment in the Typic Ustorthent. Rainfall amounts had no significant correlation 
with C mineralisation. 


Keywords: C mineralisation, C sequestration, Entisols, rainfall, temperature, SOC 
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Introduction 

Soil organic carbon (SOC) plays an essential role in the stability and fertility of soils 
(Trivedi et al., 2018) and exists in inorganic and organic forms. Soil C plays a crucial role in 
the terrestrial carbon cycle, releases nutrients for plant growth, promotes soil structure, 
biological and physical health, and is a buffer against harmful substances. These make 
maintaining optimum amounts of OC in soil paramount for sustainable agricultural 
production. 

Land management options present opportunities for good agricultural soil resource 
use, sustainable food production and sustenance of livelihoods. In this regard, optimum soil 
management under crop production is desirable for good yield and minimised adverse 
environmental effects. Previous studies have shown that land management is closely linked to 
SOC turnover (Jamala and Oke, 2013; Xiao et al., 2017), which has a far-reaching impact on 
terrestrial carbon cycles and connected ecological implications such as climate change. 

In an earlier study, Vanden et al. (2003) showed that SOC content is lost principally in 
the form of COs, following the conversion of natural land states for agricultural use. 
However, carbon-enhancing land management helps improve the carbon sink potential of 
soils (Ding et al., 2014). , the C sink or source potential of soil thus depends on soil physical 
and chemical factors influenced by land management, which impact soil functions and affect 
the activities of soil microbial populations that decompose SOM and make the nutrients 
available to plants (Guo et al., 2017). 

Climatic variables are another critical set of factors determining C turnover in soil. 
Research has shown that soil temperature and moisture affect microbial activity, surface 
runoff, land surface energy dynamics, root zone productivity, and biomass yield (Hacki et al., 
2005), all of which bear C turnover in soil. Soil water also enhance dissolved organic carbon 
(DOC) leaching from the forest floor, a vital C pool. Temperature and precipitation are 
dominant factors affecting soil organic carbon storage (Herold et al., 2014). 

Understanding the interaction between climatic factors and management options is thus 
important to design tailored land management approaches that guarantee increased food 
production and SOC buildup. This is especially important under maize production due to its 
centrality as one of the world's most important food/ feed and industrial use crops. Therefore, 
this study is designed to assess the effect of variable nitrogen, phosphorus and potassium 
compound fertilisers application rates and biochar amendment under conventional tillage 
practices for their C sequestration potential. 

Low nutrient content and accelerated mineralisation of SOM remain two significant 
constraints facing sustainable agriculture, as earlier identified by Renner (2007). In the face 
of high human population growth, especially in sub-Saharan Africa, intensified agricultural 
management has become inevitable. This reality makes studies aimed at optimising land 
management under prevailing climatic conditions pertinent. We hypothesise that biochar 
addition will lower SOC mineralisation, and N fertilisation will increase C mineralisation in 
the test soils. 

Materials and Methods 

Description of the study area 

This study was carried out within the University of Ilorin Teaching and Research 
Farm in Ilorin, Kwara State, in Southern Guinea Savannah (SGS) agroecological zone of 
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Nigeria (Figure 1). Geographically, the state is located between latitudes 7° 45' and 9° 30' 
North of the Equator and Longitudes 2° 30' and 6° 25' East of the Prime meridian. The SGS 
zone has a mean annual rainfall of 100-150 cm and a wet season lasting 6-8 months between 
March and October (Table 1). 

The climate is tropical and is characterised by double rainfall maxima tropical wet and 
dry weather. Temperature is uniformly high and ranges between 25 °C and 30 °C in the wet 
season, and dry season temperature ranges between 33 °C and 34 °C. Relative humidity in the 
wet season is between 75 to 80 %, while about 65 % in the dry season. 

The Study Area, Showing Experimental Sites 


4°39'30°E 4°40'0"E 


Map of Nigeria 
(highlighting Kwara State) 


8°28'0"N 


8°27'30"N 
8°27'30"N 


N 


. Sources: Esri, HERE, Garmin, Intermap,-increment P Corp., 
(a Experimental Area GEBCO, USGS, FAO, NPS, NRCAN, GeoBase, IGN, Kadaster 
NL, Ordnance Survey, Esri Japan, METI, Esri China (Hong 
0 125 250 500 Meters Kong), (c) OpenStreetMap contributors, and the GIS User 
Community 


4°39'30"E 4°40'0"E 
Fig. 1: Map showing the study area and experimental locations 


The state is considered one of the essential rainfed agricultural areas in semi-arid 
regions of the country and could be defined as dependent mainly on rainfall for crop 
production. Food crops produced are mostly maise, rice, sorghum, yam, cassava, water yam 
and sweet potato, which constitute the main staple foods (Ajadi et al., 2011). To a lower 
extent, the soils are formed majorly from basement complex rocks (metamorphic and igneous 
rocks) and sedimentary rocks. The metamorphic rocks include biotite gneiss, banded gneiss, 
quartzite augite gneiss and granitic gneiss. The intrusive rock includes pegmatite and vein 
quartz (Olaniyan and Ogunkunle, 2007). 
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meteorological station) 


Table 1. 1989 — 2014 monthly rainfall data for Ilorin (Source: Ilorin airport 


Years Jan Feb Mar _ Apr May Jun Jul Aug Sept Oct Nov Dec Total 
1989 0.0 0.0 56.2 96.5 106.8 239.9 111.3 1999 233.8 188.9 0.0 0.0 1233.4 
1990 0.0 8.5 0.0 158.2 130.7 129.7 150.0 1168 225.1 98.7 17.4 0.0 1035.1 
1991 17.2 17.6 37.8 67.5 340.8 1359 9.0 138.8 124.8 127.0 0.0 0.0 1326.7 
1992 0.0 0.0 2.0 46.7 166.6 90.1 134.7 75.9 129.4 193.2 0.0 0.0 838.6 

1993 0.0 102.5 79.0 35.2 176.1 2206 1004 333.6 362.8 849 0.0 0.0 1495.0 
1994 70.1 0.0 7.9 161.0 168.0 168.1 283.1 83.4 201.0 3179 289.9 0.0 1750.4 
1995 0.0 0.0 79.0 79.0 234.0 1865 140.5 245.1 272.5 133.0 11.5 12.0 1393.1 
1996 0.0 0.0 56.5 67.1 165.3 1542 96.4 1546 1998 96.3 0.0 0.0 990.2 

1997 0.0 0.0 93.1 2064 2074 375.2 1588 1076 2244 115.3 15.0 0.0 1503.2 
1998 0.0 0.8 17.1 73.9 159.6 1323 58.5 170.9 224.8 140.0 0.0 0.0 977.9 

1999 0.0 32.5 54.9 67.1 143.3. 2449 1956 85.7 265.6 182.5 22.8 0.0 1294.9 
2000 0.0 2.2 24 22.8 68.2 2622 93.4 163.9 2680 9.0 0.0 0.0 $92.1 

2001 0.0 0.0 18.0 45.0 139.0 1218 1387 44.5 176.4 60.8 0.0 0.0 744.2 

2002 0.0 0.0 59.4 163.2 572 978 180.8 182.7 144.5 176.7 6.5 0.0 1068.8 
2003 0.0 14.0 12.4 93.9 1245 360.7 123.2 130.9 1762 133.4 46.7 0.0 1215.9 
2004 18.2 33.1 4.0 67.0 260.5 1592 2114 1454 2435 981 31.6 0.0 1272.0 
2005 0.0 5.5 25.5 75.7 187.6 171.0 130.2 93.6 282.6 109.8 10.5 0.0 1092.0 
2006 1.2 16.1 27.5 106.6 163.7 2596 224.1 88.2 276.2 190.0 0.0 0.0 1353.2 
2007 0.0 7.0 22.0 115.0 306.2 227.7 205.3 189.9 236.2 162.5 0.0 2.5 1474.3 
2008 0.0 0.0 9.5 123.5 35.2 1841 394.1 329.5 3436 88.6 3.4 5.2 1516.7 
2009 6.0 0.5 19.6 219.2 112.1 1706 2495 226.0 231.0 112.0 474 0.0 1393.9 
2010 0.0 6.5 60.9 58.2 112.2 45.1 139.1 123.7 2283 2264 2.0 0.0 1002.4 
2011 0.0 19.1 48 13.5 169.1 231.8 2202 2945 334.6 163.2 0.0 0.0 1450.8 
2012 26.0 0.0 0.6 173.3 205.7 1628 2026 168.1 2008 844 0.0 7.2 1231.5 
2013 0.5 39.0 39.8 181.8 81.8 1329 107.3 17.7 202.5 154.3 0.0 114 969.0 

2014 6.3 34.2 71.0 321.4 163.8 154.7 82.1 94.9 3916 2594 0.0 0.0 1579.4 
MEAN _5.6 13.0 33.1 109.2 161.0 1854 159.2 1541 2385 1426 19.4 1.5 1234.4 


Soil Survey 

A free soil survey identified two soil types within the study area. Profile pits were dug and 
described, following United States Department of Agriculture (USDA), Soil Survey Staff 
guidelines (2017). Soil classification was done following USDA Soil Survey Staff (2014) 
guidelines. Identified soils were labelled UTRF 1 and UTRF 2 at approximately 
8°27'48.91"N 4°39'37.01"E and 8°27'21.76"N 4°39'51.53"E respectively. The soils were 
classified as Anthrodensic Ustorthent and Typic Ustorthent, respectively, based on observed 
morphological and measured physical and chemical properties. 


Laboratory Analysis 

Top (O — 130 cm) and profile horizons soil samples were air-dried, 2 mm sieved, and 
subjected to physical and chemical analyses following procedures outlined in Table 2. The 
biochar was made out of hardwood in a top-lit updraft kiln at ~350 °C. Biochar samples were 
broken, ground (using a ceramic mortar and pestle) and passed through a 0.5 mm sieve mesh 
before analyses. Analyses procedures are presented in Table 2. The soils' morphological, 
physical and chemical properties are shown in Tables 3 and 4, respectively, and biochar 
properties in Table 5. 
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Table 2. Methods of soil and biochar samples analyses 


Method of Analyses 

Parameter Soil Biochar 

Soil particle size fractions Gee and Or, 2002 

pH Thomas (1996) In 0.1 N KCl 
solution (1:10 
wt/wt) ratio 

Total Nitrogen Bremner, 1996 Enders and 
Lehmann (2012) 

Organic carbon Walkley and Black, 1934 Enders and 
Lehmann (2012) 

Available phosphorus (Bray 1-P): Bray and Kurtz (1945) Enders and 
Lehmann (2012) 

Exchangeable bases (Ca, Mg, K, 1 N ammonium acetate (NH40Ac) Enders and 

and Na) extraction Lehmann (2012) 

ECEC Juo et al. (1976) 

Exchangeable acidity Peech et al., 1962 

Percentage base saturation Summation of total exchangeable 


Moisture content 
Volatile matter 
Ash 


Fixed carbon 


bases expressed as a percentage of 
CEC 


ASTM E871-82 
ASTM E871-82 
ASTM E871-82 
ASTM E871-82 
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Table 5. Properties of biochar material 


pH (H20) 10.22 
pH (KCl) 9.24 
Electrical Conductivity (uS/cm) 5.08 
Available Phosphorus (mg/kg) 509 
Organic Carbon (%) 34.08 
Nitrogen (%) 0.28 
Ca (cmol/kg) 1.99 
Na (cmol/kg) 0.30 
Mg (cmol/kg) 1.55 
K (cmol/kg) 2.10 
ECEC (cmol/kg) 6.60 
Ash content (%) 24.3 
Moisture content (%) 4.21 
Fixed carbon (%) 48.96 
Volatile matter (%) 22.53 


Field Experiment 

The field experiment was conducted in the rainy season of 2019 (August — Nov.). It 
was laid out in a Complete Randomized Design (CRD) on each soil. Four treatments were 
used; control (C), NPK 20:10:10 at 600 kg/ha (NPK), biochar at 3.5 t/ha + NPK 20:10:10 at 
600 kg/ha (biochar+NPK), and biochar at 3.5 t/ha (B). Treatments were replicated thrice on 
plots measuring 2 x 4m (8 m’) and separated by a buffer of 1 m between plots. In NPK plots, 
200 g was applied and 400 g at 6 weeks after planting. Biochar was applied after tillage, 
before planting, and control treatments were not treated with any amendment. Quality protein 
maise (QPM), a hybrid variety, was planted immediately after amendments. Seeds were sown 
at 2 seeds per hole, 25 cm within the row and 75 cm between row spacing. Hand weeding was 
done at 21 days, 6 and 9 weeks after planting (WAP). The second NPK split application was 
made at 6 (WAP) 

Carbon Mineralisation Measurement 

Soil respiration was measured in situ using the alkali absorption method, using the 
non-flow-through steady-state (static chamber) system (Hopkins, 2006). Field chambers were 
constructed using a PVC non-reactive material following the USDA Agricultural Research 
Service guidelines in Parkin and Venterea (2010). Permanent chambers (10.6 cm diameter 
Polyvinyl Chloride (PVC) pipes of 30 cm height were installed to 20 cm depth) in the centre 
of each treatment plot at spots cleared of all organic matter residue. Three chambers were 
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used as blanks. Carbon dioxide emission was trapped in 60 mL plastic vials containing 40 mL 
0.25 M NaOH. Five millilitres (5 mL) aliquots of the trapping solution (0.25 M NaOH) were 
back-titrated to a colourless end-point with 0.15 M HCl after precipitation of carbonates by 
adding 8 mL 3M BaCl, using phenolphthalein as an indicator (Haber, 1958). CO2 was 
sampled on days 1, 3, 6 and every three days intervals till the maturity of the test crop. 
CO,-C was evaluated using the equation: 

CO2—C = (B-V)xMxE)/(AxT) 

Where: 

CO .- C = C efflux (mg m? d'); B = volume of HCI solution used for the titration of 
treatment samples (mL); V = volume of HCI solution used for the titration of blank samples 
(mL); 

M = molarity (mol L") of the HCI solution; E = carbon gram-equivalent (6 g); 

A = exposed soil surface area (m°); and T =sampling time (d). 

Meteorological Data 

Rainfall (total rainfall between gas sampling times) and soil temperature (mean 0-20 
cm) data was obtained from the University of Ilorin Teaching and Research Farm 
Meteorological station. The meteorological station is 1.6 and 0.4 km away from UTRF 1 and 
2 respectively. 

Data Analysis 

Data collected on C mineralisation was subjected to analysis of variance (ANOVA) at 
p = 0.05. In contrast, rainfall, soil temperature and carbon mineralisation data were subjected 
to Pearson Correlation (r) analysis to explore the relationships between the variables. 
Calculations were done in SPSS 20" Edition. 

Results and Discussion 

Variation in C Mineralization and CO, Emission Treatments and Soils 
Table 6 shows the results of the analyses of variance of means of experimental factors. The 
highest CO2 emissions were recorded under NPK treatment (79.856 mg), while the least was 
in the biochar alone treatment at 74.845 mg. However, the differences in the values for all 
treatments were not significant at p<0.05. CO2 emission rate followed a similar trend with that 
of total emissions, having no significant difference among treatment applications. Soils, 
however, exerted significant variations both in total CO 2 emission and emission rate. The 
highest values were recorded in UTRF 2, with observations of 82.351 mg and 1.183 mg m?d 
' for total CO: emission and mean emission rates, respectively. These values were 
significantly higher than those recorded in UTRF | at p<0.01 and p<0.05, respectively. No 
significant interactions were recorded between experimental factors. 

The observations recorded showed that sole NPK treatment and that incorporated with 
biochar resulted in higher C mineralisation. This observation is different from the findings of 
Olaniyan et al. (2020) in their evaluation of biochar C stability in soils around the same 
geographical location as in this study. In their research, they found that SOC and biochar C 
mineralisation was lower in soils of higher organic N, and they attributed their finding to the 
possible toxicity effect of high N on soil C degradation microbial communities. However, the 
N contents of the soils in the study are much higher than those in this research. This may 
explain the variation observed. At lower soil N levels, N addition may increase microbial 
population and activities rather than inhibit it, which will result in higher N mineralisation. 
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The finding also differs from that of Perveen et al. (2019) but is consistent with that of 
Wang et al. (2015), which reported that N enrichment resulted in higher CO2 release in 
croplands. The variation in C mineralisation responses between soils can be attributed to 
differences in properties. While both soils are comparable in properties, especially soil 
physical (particle size fractions) with a similar range in sand silt and clay components, UTRF 
1 had a higher OC amount. This can be attributed to the propensity of low OC soils to better 
adsorb soil C on soil clay mineral surfaces, which makes it less free for microbial degradation. 
Kimetu et al. (2009) also reported that amended C is less accessible for microbial utilisation 
in low N and C soils. 

Table 6. Variation in means of CO2 emissions in treatments and soils 


Total CO, emission Mean CO, emission rate 


(mg) (mg m” d’') 
Treatment 
B 74.845 1.072 
Biochar+NPK (BN) 78.964 1.134 
C 77.494 1.125 
NPK 79.856 LAS] 
F-value 0.512 0.696 
p-value 0.680 NS 0.568 NS 
Soil 
UTRF 1 73.228 1.061 
UTRF 2 82.351 1.183 
F-value 8.870 7.828 
p-value 0.009 ** 0.013 * 
Interaction 
TxS 0.510 NS 0.382 NS 


*: significant at p<0.05, **: significant at p<0.01, NS: not significant 


Relationship Between Soil Temperature and C Mineralization in Soils and 
Treatments 

Fig. 2 shows topsoil (0 — 20 cm) temperature and CO 2 emission rate trends in 
treatments and soils. There was no clear trend in C mineralisation and emissions in both soils 
with soil temperature. The general direction was highly fluctuating emission rates at the start 
of incubation before a stabilisation at about day 45 and 57 in UTRF | and 2, respectively. 
This observation largely conforms with previous findings of a more variable emission at the 
start of incubation and a steady progression with time (El-Naggar et al., 2015; Olaniyan et al., 
2020). 

However, correlation analyses showed that the relationship between soil temperature 
and C mineralisation differed between treatments and soils. In Table 7, which presents the 
results of correlation analysis on an individual soil basis, there was a significant positive 
relationship between soil temperature and SOC mineralisation under biochar (B) treatment 
(r=0.44; p=0.03; N=24) in UTRF 2, whereas the relationship was not significant in UTRF 1, 


Bulgarian Journal of Soil Science® 2022 Volume 7. Issue 1. www.bsss.bg 42 


Isimikalu et al., 2022 


and the correlation was negative. No other correlation for the different treatments was 
substantial. 

A general correlation analysis (Table 8) irrespective of soils also showed that soil 
temperature had a significant positive correlation with C mineralisation rate under biochar 
treatment (r=0.29; p=0.04; N=24), and the same observation as in the soil-wise analysis was 
observed, with no significant correlation under other treatments. Although insignificant, the 
relationship between soil temperature and C mineralisation was positive under all 
experimental treatments. This finding aligns with that of Benbi and Khosa (2014), in who 
they found that incubation medium temperature significantly influenced C mineralisation. 
Although they also noted that the chemical composition of the C source also determines how 
much effect temperature has on C mineralisation. This also applies in the current study. The 
variation in inherent C sources in the soils can be attributed to the difference in the 
relationships observed between soil temperature and C mineralisation in the two soils. 


Treatments CO, Emission Rates and Soil Temperature Trend in UTRF 1 
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Fig. 2. Soil temperature and CO2 emission rate trends in soils and treatments 
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Table 7. CO2 emission rate, rainfall and soil temperature relationships in soils and 


treatments 
Soil Treatment Rainfall (mm) Top (20 cm) Soil Temp 
Emission Rate (°C) 
(mg m” da’) 
UTRF NPK r -.080 .040 
1 P eae 853 
c r 144 059 
p 503 .784 
BN r 110 118 
P .610 582 
B r -.100 -.017 
P 642 939 
UTRF NPK r .170 .295 
2 P 427 161 
C r -.055 .287 
P 197 174 
BN r 32d 238 
P 123 .263 
B r 211 444° 
P Be Se .030 


N=24, *. Correlation is significant at the 0.05 level (2-tailed), CN: control treatment; B: biochar treatment; 
NPK: NPK treatment; BN: biochar-NPK treatment 


Table 8. General CO 2 emission rate, rainfall and soil temperature relationships in 
experimental treatments 


Treatment Top (20 cm) Soil Temp Rainfall (mm) 
(°C) 
B Emm rate (mg m? d') R 0.289* 0.096 
P 0.046 0.516 
BN Emm rate (mgm*d')  R 0.264 0.224 
P 0.070 0.127 
C Emm rate (mg m? d') R 0.216 0.058 
P 0.140 0.698 
NPK Emm rate (mg m’d') R 0.261 0.081 
P 0.074 0.586 


N=24, *. Correlation is significant at the 0.05 level (2-tailed), CN: control treatment; B: 
biochar treatment; NPK: NPK treatment; BN: biochar-NPK treatment 
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Relationship Between Soil Moisture and C Mineralization in Soils and 
Treatments 


Fig. 3 shows the trends of soil moisture and C mineralisation in this study. As with soil 
temperature, there was no strong consonance between soil moisture and C mineralisation 
fluctuations. This is reinforced by the correlation analysis carried out. The results showed that 
rainfall amount had no significant correlation with C mineralisation and CO, emissions under 
treatments in both soils. Also, the direction of relationships differed slightly between 
treatments, with NPK and B treatments having a weak negative correlation in UTRF | and the 
control (C) treatment in UTRF 2. 
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Fig. 3. Rainfall trend and treatments CO2 emission rates in soils and treatments 


This finding supports models that have shown that C mineralisation declines under 
elevated moisture. In their study, Benbi and Khosa (2014) found that the flooding of 
experimental soils reduced the mineralisation of added C. In a separate study, Huang and Hall 
(2017) reported an initial (25 days) suppression of C mineralisation, following water 
saturation above field capacity. However, after this initial period, they reported that elevated 
moisture stimulated cumulative CO2 and CH4 losses up to over 150 % compared to the 
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control. Stable C isotopes examination in their study showed that C losses were derived from 
C3-C release under Fe reduction at high soil moistures. 

The observations in this study suggest that soil C losses with soil water may be more 
connected to the leaching DOC. Studies have shown that moisture-induced high C 
mineralisation results in moisture levels that optimise soil oxygen supply and carbon substrate 
diffusion (Linn and Doran, 1984). As this study was carried out during peak rainfall in the 
study area, which results in high water saturation above field capacity, the saturating may 
have resulted in the C mineralisation observations recorded. 


Conclusion and Recommendations 

In examining the influence of soil management, temperature and moisture on C 
mineralisation in this study, the finding showed that assessed management options did not 
significantly affect C mineralisation. However, the two soils evaluated differed in C 
mineralisation responses. This observation could be attributed to a variation in soil properties, 
most likely OC, since the soils were mainly similar in other properties except in OC gradient. 
Findings showed that C mineralisation is more sensitive to soil temperature under sole biochar 
treatment. Although this observation did not cut across both soils studied individually, 
analyses showed that higher temperatures increased C mineralisation under exclusive biochar 
treatment. The outcome of this study can be applied in timing biochar application to the soil to 
maximise its carbon sequestration potential. 
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Abstract 

Crude oil contamination of farmlands and rivers has been a major concern in Nigeria, 
especially in places like Ogoni land in Rivers State and some other oil producing states. It has 
severely impacted on the economic and health status of people in that region. It is important to 
find ways to ameliorate the harmful effects of crude oil contamination on both terrestrial and 
aquatic lives. In this study, we investigated the bioremediation profile of sandy and silty clay 
soils contaminated with crude oil. Contaminated soils were treated with inorganic fertilizer 
(NPK 15:15:15) and urea to improve the density of soil microbes and enzymes. The total 
petroleum hydrocarbon (TPH), total heterotrophic bacteria (THB), total hydrocarbon utilizing 
bacteria (THUB), total heterotrophic fungi (THF), total hydrocarbon utilizing fungi (THUF), 
and the activities of soil enzymes such as laccase, peroxidase, lipase, and catalase were 
investigated. The results indicated that treatment of crude oil contaminated soils with NPK 
and urea increased the activities of soil enzymes and enhanced total soil microbial load that 
play active roles in biodegradation of petroleum products, thereby suggesting that NPK and 
urea treatments of soils have the potential of decreasing crude oil contamination and 
remediating contaminated soils. Furthermore, this study showed that NPK and urea decreased 
the total petroleum hydrocarbon (TPH) values of both sandy (92.5%) and silty clay (80.5%) 
soils. The findings from this study, suggest that NPK and urea application to contaminated 
soils could serve as beneficial remediation agents for quick mitigation of the hazardous effects 
of crude oil contamination on farmlands. 


Keywords: Crude oil contamination, (Bio)remediation, Soil microbes and enzymes, 
NPK. 


Introduction 

Nigeria is one of the countries with large stores of crude oil, especially in the Southern 
part, and due to continuous oil drilling and exploration, there are many sites, including 
residential areas, farmlands, rivers, and the surrounding atmosphere that have been polluted. 
The pollution has affected both human and animal lives, through introduction of harmful 
petrochemical toxicants into the environment. The toxicants cause serious damage to various 
tissues and organs in the human body, terrestrial animals, and aquatic lives. Also, it affects 
farmlands and results in poor yield of agricultural products, and as a consequence lead to poor 
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economic status of the inhabitants. There are many petroleum industries that operate in the 
Southern part of Nigeria due to their rich oil store, and their operations, over the years, have 
led to the pollution of these territories and have rendered many of them uninhabitable. They 
have caused tremendous hardship to the communities and have affected the ecosystem 
adversely resulting in a loss of biodiversity (Akpofure et al., 2000). Oil spillage on farmlands 
usually reduces yields and destroys soil fertility until there is a proper clean-up of the 
farmland to restore the natural ecosystem processes and reestablish stability (Sada and 
Odemerho, 2000). 

Different methods have been employed to remediate crude oil polluted sites, reduce 
the damage caused, and restore the natural processes that had been altered by the pollutants. 
Some of the methods include natural method, also known as natural attenuation; physical 
methods such as booming, skimming, washing, tilting, and sediment relocation; chemical 
methods such as the use of detergents, dispersants, de-emulsifiers, and surface film chemicals. 
Another viable method that has been applied in cleaning up polluted farmlands is 
bioremediation, which is a biological method. It is the use of biological agents which are 
capable of utilizing the toxic petrochemical agents as energy source, and in the process 
degrade them to non-toxic forms (USEPA, 1989). 

Microorganisms and some plants have been reported as key players in the 
bioremediation processes. When plants are used to bioremediate a polluted soil, it is termed 
phytoremediation. Microorganisms such as bacteria and fungi that live ubiquitously and 
utilize large number of petrochemical compounds as energy source are essential agents in the 
bioremediation processes. They can degrade the contaminants and use them as source of 
energy because they are supplied with various enzymes that catalyze the breakdown of these 
compounds. Fungi produce extracellular enzymes directly unto the contaminants and degrade 
them to usable and less toxic forms which are then absorbed into the cells. In the case of 
bacteria, the contaminants are first absorbed into the cells before they are acted upon by the 
degrading enzymes. The combined action of fungi and bacteria results in a clean-up of the 
petrochemical toxicants and bioremediation of the polluted soils. Some of the degrading 
microorganisms that have been identified include Aspergillus niger, Pseudomonas, and 
Penicillium (Cooney, 1994). The supply of required nutrients to these microorganisms may 
enhance their bioremediation capacity of polluted soils. Therefore, in this study, crude oil- 
polluted soils were treated with NPK 15:15:15 and urea fertilizer as a way of supplying the 
resident microorganisms with nutrients to enhance their bioremediation capacity. Hence, the 
use of nutrients such as NPK and urea fertilizers on bioremediation of crude soils in the 
Niger-Delta region of Nigeria was investigated. 


Materials and Methods 

2.1. Collection of Soil Sample 

Soil samples were collected from an unpolluted site in Nigeria. The soil samples were 
collected from a depth of 0-50 cm in a sterile polythene bag. 


2.2. Determination of Physicochemical Properties 


Two soil types (Sandy and Silt clay) were used for this study, and their physical and 
chemical properties were determined using reagents of analytical grade. The following were 
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the physicochemical properties determined: Available phosphorus, Total nitrogen, Total 
Organic Carbon (TOC), pH, bulk density, porosity, water holding capacity and electrical 
conductivity (EC). These properties were determined three days (0, 18, and 36 day). 

2.3. Soil Temperature Measurement 

A pilot hole was made in the soil using a screw driver, and the thermometer was 
inserted into the hole making sure it’s firm and kept for several minutes until a stable 
temperature was achieved (Hewins et al., 2016). The temperature of the soil was taken at the 
site. 

2.4. Soil pH Determination 

The method of McLean (1982) was adopted. The soil was air-dried and 10 g was 
placed in a 25 mL beaker. Deionized water (10 mL) was then added to the soil sample and 
thoroughly mixed for 5 sec with glass rod, and the suspension was allowed to stand for 30 
min. The pH was measured by inserting the electrode of calibrated pH meter into soil 
suspension. 

2.5. Determination of Bulk Density 

An empty universal bottle was weighed and filled with oven-dried soil sample up to 
the brim by tapping, and the weight was recorded. The volume of the bottle was determined 
using a burette (FAO, 1980), and the bulk density was calculated as follows: 


: : weight oven dry soil 
Soil Bulk Density (g/cm?) ———r——rro 
volume of soil 


2.6. Determination of soil moisture content 

Gravimetric method was used for this analysis (Reynolds, 1970). It involved the 
determination of mass difference between the wet soil and soil oven-dried at 105°C. An empty 
crucible was weighed and the mass was recorded as M, (g). the wet soil sample was 
positioned in a crucible and weighed immediately and the mass was recorded as Mp (g). The 
crucible containing the wet sample was placed in an oven and dried at 105°C until a constant 
weight was achieved. The sample was removed from the oven and placed in a desiccator to 
cool. The crucible with the oven-dried soil was weighed and the mass was recorded was M3 
(g). The formular below was used to calculate moisture content: 


M2-M3 
x100 
M3-M1 


Where M, = Mass of empty crucible (g), Mz = Mass of crucible and wet soil (g), M3 = 
Mass of crucible and oven dried soil (g) 

2.7. Water Holding Capacity of Soil 

A small cup with perforated bottom was weighed and recorded. Filter paper was 
placed at the foot of the cup and the weight was taken. Soil sample (70 g) was then weighed 
into the cup on a retort stand and 250 mL measuring cylinder with funnel was inserted into it. 
Then, 50 mL of deionized water was added and it was allowed to drain overnight. The drained 
water was read and recorded (FAO, 1980). Water holding capacity of soil was calculated as 
follows: 


Moisture Content (%) = 


Volume of water retained 
WHC = ————————_ * 100 


Volume of sample 
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2.8. Soil Porosity 

Air-dried soil was added to 500 mL measuring cylinder and 100 mL of deionized 
water was carefully and slowly added to the soil sample inside the measuring cylinder until it 
reached the top of the soil. The volume of water used was recorded (FAO, 1980), and porosity 


was calculated as follows: 


: Amount of water added to soil sample 
Porosity = ————_______—_——— * 100 
Total Soil Sample 


(Porosity (%) = 100- (bulk density/particle density) x 100) 


2.9. Determination of Total Nitrogen in Soil 

Soil sample (2 g) was weighed into an 800 mL Kjeldahl digestion flask, deionized 
water of 10 mL was added, followed by addition of 5 g of Kjeldahl catalyst mixture and 15 
mL of 18.6 M sulphuric acid. The flask was cautiously heated on digestion stand until 
frothing stopped. The heat was increased until digest was clear, the flask was cooled and 
deionized water was slowly added and shaken. Three chips of granulated zinc were added and 
30 mL of 2% boric acid were measured into 250 mL conical flask and placed under a 
condenser. Then, 75 mL of 10 M NaOH was slowly added into the digest in Kjeldahl flask 
and immediately connected to distillation apparatus. The distillate with NH3 liberated was 
titrated with a standard acid along with four drops of mixed indicator (FAO, 1980). Total 
nitrogen was calculated as follows: 
Total Nitrogen (%N) = 
(Equivalents of acid added to sample — equivalents of acid added to blank)(14.01)(100) 

Sample weight (g) 


2.10. Available Phosphorus in Soil 

Air-dried soil (2.5 g) was placed in a clean and dried 125 mL polyethylene bottle, 
followed by the addition of 50 mL of 0.5 M sodium bicarbonate (pH 8.0) with 
polyacrylamide, which was used as the extracting solution from a dispenser. The bottle was 
placed in a reciprocating shaker and shook for 30 min and filtered into test tube with filter 
paper. Technicon II Autoanalyzer was used to determine the available phosphorus in clear 


filtrate by colorimeter. Soil phosphorus was calculated as: 
AXBxCxM 
E 
Where, A = Sample extract reading (mg I-1), B = Extract volume (ml), C = Dilution, if 


performed, M = Moisture correction factor, E = Sample weight (g) 


Soil Phosphorus (mg/kg) = 


2.11. Total Petroleum Hydrocarbon (TPH) 

Sandy and silty clay soil samples (5 g) contaminated with 5000 ppm of crude oil or co- 
treated with NPK and urea at CPN ratio of 100:2:0:2 were extracted three times with 15 mL 
of hexane each. The three extracts (three) were pooled and dried by solvent evaporator at 
room temperature under gentle nitrogen stream in a fume hood. Residues were weighed along 
with solvent blank beaker. The difference of the weight was calculated as the total petroleum 
hydrocarbon content (Chandra et al., 2013). 

2.12. Gas Chromatographic Analysis of Soil 
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Five grams (5 g) of the soil sample was mixed with 5 g of anhydrous sodium sulphate in a 25 
mL vial, 15 mL of n-hexane was added and vortexed for 10 min. The resulting soil suspension 
was filtered using a 0.45 um Teflon filter. The filtrate was analyzed using a gas 
chromatography with flame ionization detector (GC-FID) (HP 7890). The injector and 
detector temperatures were programmed at 250°C and 350°C, respectively. The initial oven 
temperature was maintained at 50°C, ramped at a rate of 5 - 280°C per min, and held for 6 
min. The injection volume for both sample and standards were | pL. 

2.13. Soil Enzymes Assays 

2.13.1. Laccase Activity 

Laccase activity was determined using pyrogallol as substrate (Allison and Jastrow, 
2006). Briefly, 1 g of soil sample was weighed into a 250 mL conical flask and 50 mL acetate 
buffer (50 mM, pH 5.0) was added. The flask was incubated at room temperature for 1 h with 
vigorous shaking every 20 min. The volume of the buffer was increased to 125 mL and the 
flask shaken vigorously. Aliquot of 10 mL of the soil suspension was transferred to a 
centrifuge tube and centrifuged at 4000 rpm for 10 min. The supernatant obtained was used 
for the enzyme assay. The test experiment contained 2 mL of the supernatant (soil suspension) 
and | mL substrate (25 mM pyrogallol). This was incubated in the dark at room temperature 
for 1 h. A sample control containing 2 mL supernatant and | mL buffer and also substrate 
control containing 1 mL substrate and 2 mL buffer were treated as the test experiment. The 
absorbance was measured at 460 nm using UV — Visible spectrophotometer. The buffer 


solution was used as blank. Laccase activity was calculated as follows: 
AXxV1 
E xV2 x T XW 
Where, A = Net absorbance = Test — sample control — substrate control, V; = Volume of 


Laccase activity (umol/h/g) = 


buffer used, E = Molar extinction coefficient for pyrogallol = 4.2/umol, V2 = Volume of soil 
suspension, T = Substrate incubation time, W = Mass of soil sample 

2.13.2. Peroxidase Activity 

Peroxidase activity was estimated using pyrogallol as substrate (Allison and Jastrow, 
2006). The procedure described for laccase activity was followed, but with the addition of 0.2 
mL of 0.3% hydrogen peroxide to the test, sample control and substrate control experiments 
respectively. 

2.13.3. Lipase Activity 

Lipase activity was determined using the method described by Parry et al (1966). 
Briefly, 1 g of soil sample was weighed into a 250 mL conical flask and 100 mL phosphate 
buffer (50 mM, pH 7.4) was added. The conical flask was shaken vigorously and then 
incubated at room temperature for 1 h with intermittent shaking. Thereafter, 10 mL of the soil 
suspension was transferred to a centrifuge tube and centrifuged at 4000 rpm for 10 min. The 
supernatant was used for the lipase activity assay. The test experiment was made up of 1 mL 
supernatant and 5 mL 10% olive oil and gum Arabic emulsion in a conical flask. The flask 
was incubated for 1 h on a shaker with the speed adjusted to high and temperature set at 40°C. 
At the end of the incubation, 1 mL 95% ethanol was added to the conical flask and shaken to 
stop the reaction. The free fatty acid released was determined using a UV spectrophotometer. 
Control experiment was also one with 1 mL buffer solution. 
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Lipase activity was calculated as follows: 
(VS—VB)xN x1000 
S 
Where Vg = Volume of NaOH used against control flask, Vs = volume of NaOH used against 


experimental flask, N = Normality of NaOH, S = Volume of enzyme extract, 


Lipase activity (U/g) = 


2.13.4. Catalase Activity 

Catalase activity was determined according to the method of Cohen et al (1970). 
Briefly, soil sample supernatant was obtained as described for lipase activity. The supernatant 
(1 mL) was transferred to a test tube and 5 mL of 0.3% hydrogen peroxide was added. The 
content of the test tube was mixed by shaking and was then incubated at room temperature for 
5 min. The reaction was stopped by the addition of 1 mL of 6 N sulphuric acid and 0.01 M 
KMn0O, was added within 3 sec. The absorbance was measured using a spectrophotometer at 
480 nm within 30-60 sec. Control experiment was performed using deionized water. Catalase 
activity was calculated as follows: 
(Absorbance) xv x1000 

M xv XW 

Where, Absorbance = absorbance (control) — absorbance (test), V = total volume of reaction 


Catalase activity (U/g) = 


mixture, M = molar extinction coefficient = 40.0, v = volume of supernatant used, W = mass 
of soil sample, t = reaction incubation time = 5 min. 


2.14. Effect of Abiotic Treatment on Crude Oil Biodegradation in Soil 
The ex-situ remediation method of Saterbak et al (2000) was adopted with modification. 
Briefly, each soil sample (3 kg) contaminated with crude oil (3000, 5000, and 8000 mg/kg) 
was weighed into 7 plastic buckets for each soil type (sandy soil and silty clay). Moisture 
content (MC) was adjusted to 60-80% of water holding capacity (WHC) for all group. Each 
container was agitated every three days and same with moisture adjustment. The experiment 
lasted for 36 days with the following determinations done every 6 days: Microbial growth 
were analyzed using, THB, THUB, THF, and THUF as markers. 

2.15. Statistical Analysis 
Experimental data were analyzed using IBM SPSS statistics 23 software for Windows. All 
data are expressed as mean + SEM (standard error of mean). ANOVA was used in comparing 
the mean followed by Duncan’s Multiple Range (DMRT), Post Hoc test. Student’s t test was 
used to compare means when only two means were involved. Statistical significance was 
taken as p < 0.05. Correlation analysis was done to establish relationship between TPH and 
the other parameters evaluated (THUF, THF, THB, THUB, laccase, peroxidase, catalase, 
lipase activities) respectively. 

Results 

3.1. Physicochemical Properties of Silty Clay and Sandy Soils Contaminated with 
Crude Oil 

The results below show the physical and chemical properties of crude oil- 
contaminated sandy soil and silt clay. 
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Table 1. Physicochemical data on silty clay and sandy soils contaminated with crude 
oil. Values are expressed as mean + SEM (n = 3, at 0, 18, and 36 day). 


Silty clay contaminated 


Sandy soil contaminated 


Parameters “Control 3000-5000 ~-~« 8000 ~—~*| Control 3000 ~~—«5000+~—~—~8000 
ppm ppm ppm ppm ppm ppm 
pH 6.84 + 5.70+ 5.50+ 5.70+ 6.33 + 6.7140 6.514 6.80+ 
0.05 0.01 0.04 0.03 0.31 .13 0.36 0.63 
Total nitrogen 0.49 + 0.18+ 0.22+ 0.27+ 0.43 + 0.1740 0.214 0.22+ 
(mg/kg) 0.00 0.03 0.01 0.49 0.76 51 0.11 0.11 
Available 7.10 + 4.42+ 4.48+ 4.53+ 6.10 + 3.42+1 3.48+ 3.53+ 
phosphorus 2.38 1.10 1.15 2.20 2.09 25 1.32 1.10 
(mg/kg) 
TOC (%) 0.07+0.4 0.104 0.21+ 0.43+ 0.06 + 0.1040 0.29+ 0.43+ 
+0.2 0.21 0.12 0.16 12 0.01 0.19 
Potassium 8.91 + 6.55+ 6.59+ 6.61+ 6.93 + 4.5441 457+ 4.59+ 
2.01 1.98 1.52 1.31 1.65 11 1.21 1.41 
TPH (mg/kg) <1 7568 7601 7608 <1 7568 7601 7608 
Moisture 0.30+ 0.32+ 0.74+ 1.32+ 0.25+ 0.2840 0.69+ 1.27+ 
content (%) 0.02 0.04 0.01 0.05 0.08 .O1 0.03 0.06 
Soil bulk 1.09+ 2.87+ 2.92+ 2.02+ 1.69+ 2.8140 3.944 4.08+ 
density (g/em*) 0.02 0.12 0.15 0.98 0.76 23 0.31 0.68 
Conductivity 180+ 3.10 190+ 170+ 190+ 160+ 2.20 170+ 150+ 1.89 170+ 
(us/cm) 2.87 1.99 2.55 2.01 1.77 
Porosity 0.51+ 0.36+ 0.31+ 0.29+ 0.31+ 0.3340 0.36+ 0.394 
0.04 0.07 0.04 0.07 0.02 01 0.09 0.07 
Water holding 0.19+ 0.19+ 0.19+ 0.19+ 0.16 0.1640 0.16+ 0.16+ 
capacity 0.00 0.01 0.01 0.01 + 0.03 02 0.03 0.02 
Temperature 27 
°C 
Particle size 
distribution 
* Silty clay % 72.1 
* Sand % 
93 


3.2. Physicochemical Properties of Contaminated Silty Clay and Sandy Soils 


Treated with NPK and Urea Fertilizers. 
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The results below show the physical and chemical properties of contaminated sandy 
soil and silty clay treated with NPK 15:15:15 and urea at a CNP (carbon, nitrogen, and 


phosphorus) ratio of 100:2:0:2. 


Table 2. Impact of NPK and urea fertilizers on the physicochemical properties of silty 


clay and sandy soils contaminated with crude oil. Values are expressed as mean + SEM (n = 
3, at 0, 18, and 36 day). 


Parameters Treated silty clay Treated sandy soil 

3000 ppm 5000 ppm 8000 ppm =| 3000 ppm )=5000 ppm 8000 ppm 
pH 5.69+0.31 6.18+0.33 6.44+0.10 8.3640.21 8.2840.19 8.44+0.23 
Total Nitrogen 0.31+0.01 0.37+0.04 0.42+0.02 0.28+0.11 0.3140.12 0.3940.22 
(mg/kg) 
Available 6.19+1.09 6.38+1.23 6.7741.32 5.14+0.98  5.3340.87 5.74+0.92 
Phosphorus 
(mg/kg) 
TOC % 1.26+0.21 2.02+0.12 4.48+0.31 1.26+0.23 2.02+0.22 4.48+40.31 
Potassium 72.4+4.01 73.543.81 74.843.17 70.642.12 71.442.02 72.1+1.87 
THC (mg/kg) 634+3.02 619+3.41 540+3.44 634+3.11 61942.45  540+3.65 
Moisture 3.81+£0.33 3.93+0.23 4.02+0.31 3.7640.43  3.88+0.14 3.97+40.22 
content (%) 
Soil bulk 3.17+40.12 3.28+40.22 3.39+0.12 3.722019 3.8640.21 3.9140.32 
density (g/em*) 
Conductivity 230+4.40 20043.98 22843.99 210+3.01 180+3.41  208+3.06 
(us/cm) 
Porosity 0.24+0.01 0.34+0.02 0.35+0.01 0.42+0.04 0.44+0.02 0.43+0.01 
Water Holding 1.24+0.15 1.24+0.57 1.26+0.43 1.22+0.34  1.22+0.32 1.24+0.48 
Capacity 
Temperature °C 26 P| 


Particle Size 
Distribution 
*Silty clay (%): 
72.1 
*Sandy 
(%): 93 


soil 
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3.3. Effect of NPK and Urea on Soil Enzymes Activity 

To investigate the effect of NPK and urea fertilizer treatment on enzymes activities of 
contaminated soils. Two different soil types, silty clay and sandy, contaminated with crude 
oil, were treated with NPK and urea fertilizers. The results indicated that the activities of soil 
enzymes were affected by applied treatment during the 36 days period. An increased activity 
was observed for laccase from days 0-18 on both soil types across the treated groups. 
However, laccase activity in sandy soil was higher compared to that of silty clay (Fig. 1.A). 
Peroxidase activity increased in 12 days from the start of the experiment across treated groups 
in both soils. Peroxidase activity in silty clay was higher than that of sandy soil (Fig. 1.B). 
Furthermore, an increased activity was recorded for catalase in sandy soil and silty clay with 
the highest peak at day 18. Catalase activity in sandy soil was higher than that of silty clay 
(Fig. 1.C). Also, lipase activity increased from days 0-12 in sandy soil while for silky clay, the 
highest activity was on day 24 (Fig. 1.D). Taken together, these results suggested that 
treatment of crude oil contaminated soils with NPK and urea fertilizers has the potential of 
enhancing soil enzymes activity that function in degrading crude oil contaminant. 
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Fig. 1. (A) Effect of NPK and urea fertilizers on soil enzymes activities. Effect of NPK 
and urea on soil laccase activity. Values are expressed as mean + SEM (n = 3). P < 0.05. 
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Fig. 1. (B) Effect of NPK and urea fertilizers on soil enzymes activities - Effect of 
NPK and urea on soil peroxidase activity. Values are expressed as mean + SEM (n = 3). P < 
0.05. 
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Fig. 1. (C) Effect of NPK and urea fertilizers on soil enzymes activities. Effect of NPK 
and urea on soil catalase activity. Values are expressed as mean + SEM (n = 3). P < 0.05. 
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Fig. 1. (D) Effect of NPK and urea fertilizers on soil enzymes activities. Effect of NPK 
and urea on soil lipase activity. Values are expressed as mean + SEM (n = 3). P < 0.05. 


3.4. Effect of NPK and Urea Treatment on Soil Microbial Load 


To investigate the effect of NPK and urea fertilizers on soil microbial load. Sandy and 
silty clay soils were treated with NPK 15:15:15 and urea, and the microbial load was 
monitored and reported as CFU/g during the 36 days period of treatment. An increase in 
amount of total hydrocarbon utilizing fungi (THUF) and total heterotrophic fungi (THF) were 
observed from day O and peaked on day 12 for both sandy and silty clay soils before 
decreasing thereafter in all the treatments (Fig. 2.A and B). The same pattern of growth was 
observed for total hydrocarbon utilizing fungi (THUB) and total heterotrophic fungi (THB) 
with maximum peak on day 18 (Fig. 2.C and D). These results indicated that treatment of 
crude oil contaminated soils with NPK and urea fertilizers enhanced the total soil microbial 
load. Increased soil microbial load results in increased degradation of crude oil contaminant in 
the soils and hence helps to ameliorate the toxic effects of crude oil on soils. 
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Fig. 2. (A) Effect of NPK and urea on soil microbial load of soil contaminated with 
crude oil. Effect of NPK and urea on THUF of sandy and silty clay soils. Values are 
expressed as mean + SEM (n = 3). P < 0.05. 
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Fig. 2. (B) Effect of NPK and urea on soil microbial load of soil contaminated with 
crude oil. Effect of NPK and urea on THF of sandy and silty clay soils Values are expressed 
as mean + SEM (n = 3). P < 0.05. 
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Fig. 2. (C) Effect of NPK and urea on soil microbial load of soil contaminated with 
crude oil. Effect of NPK and urea on THUB of sandy and silty clay soils. Values are 
expressed as mean + SEM (n = 3). P < 0.05. 
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Fig. 2. (D) Effect of NPK and urea on soil microbial load of soil contaminated with 
crude oil. Effect of NPK and urea on THB of sandy and silty clay soils. Values are expressed 
as mean + SEM (n = 3). P < 0.05. 
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3.5. Effect of NPK and Urea Treatment on Total Petroleum Hydrocarbon 


The effect of NPK and urea fertilizers on oil degradation at 5000 ppm on sandy soil 
and silty clay was investigated. A significant decrease in total petroleum hydrocarbon (TPH) 
concentration was observed in both soil types. Comparatively, the rate of TPH reduction in 
sandy soil was higher than that of silty clay as the former having 92.90% TPH reduction and 
the latter having 80.46% after the 36 days enhanced treatment. 
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Fig. 3. Effect of NPK and urea on TPH of soils contaminated with crude oil. Values 
are expressed as mean + SEM (n = 3). P < 0.05. 


Discussion 

Petrochemical discharge from crude oil industries on farmlands is one of the major 
problems faced by some parts of Nigeria. These petroleum products are harmful to body 
tissues and organs resulting in life-threatening disease conditions. Plants and aquatic lives are 
affected as well by the toxic nature of these substances. Therefore, there is an urgent need to 
clean up the pollution to restore normalcy of the ecosystem. This study investigated the 
additional effects of nutrients, in the form of NPK and urea fertilizers, on the bioremediation 
of crude oil-polluted soils. The results indicated that crude oil polluted soils had higher 
electrical conductivity, bulk density, and water holding capacity as compared to the 
unpolluted soil. This was in agreement with the report of other investigators (Barua et. al., 
2011). However, other properties as pH, total organic carbon, potassium, moisture content, 
and phosphorus were higher in the uncontaminated soil compared to the contaminated sample 
(Table 1 and 2). Crude oil in the soil is known to affect the properties and productivity of soil 
(Bosma et al., 1997; Michel and Fingas, 2016). The sandy nature, porosity and the pH of the 
contaminated soil supports leaching and possible contamination of the ground water (Bosma 
et al., 1997; Michel and Fingas, 2016). Inorganic nutrient sources like NPK 15:15:15 and urea 
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with nitrogen content of 46.6% was used in this study, because the quantification of required 
nutrient level in bioremediation studies is easier when inorganic nutrients are used (Suja et al., 
2014; Shahi and Aydin, 2016). 

Laccase and peroxidase are key extracellular oxidative enzymes produced by 
oleophilic microbes (Bach et al., 2013; Godoy et al., 2016). They were assayed in this study 
and a negative relationship was established between them (oxidative enzymes activity) and 
TPH using Pearson correlation. This indicated that increase in oxidative enzyme activity 
corresponded with decrease in TPH in soils (Fig. 1.A and B). However, the correlation 
between enzyme activity and microbial load was significant and positive indicating a 
corresponding increase of both enzyme activity and microbial load during the degradation 
process in both sandy soil and silt clay. The activities of catalase and lipase were used to 
monitor the reduction in toxicity level following oil degradation in soil (Margesin and 
Zimmerbauer, 1999; Riffaldi et al., 2006; Mahmoud et al., 2016), because they provide useful 
information about the soil quality. There was an increase in activity of both catalase and 
lipase in all treated variants. Catalase activity in soil decreased with increasing oil 
concentration (Fig. 1c). Wu et al (2016) corroborated this finding and affirmed the suitability 
of catalase activity for monitoring the bioremediation of crude oil. Other researchers including 
Achuba and Peretiemo-Carke (2008), Ajao et al (2011), and Achuba and Okoh (2014) have 
all reported increase in catalase activity during bioremediation as a mark of increased 
microbial activity. Similar result was observed in this study when catalase activity was 
monitored in treated sandy soil and silty clay. Lipase catalyzes the hydrolysis of glycerides at 
water-lipid interface to enable the hydrolysis products of lipids to be absorbed and used for 
metabolic activities by the microbes (Margesin and Zimmerbauer, 1999; Mahmoud et al., 
2016). It was used by Margesin and Schinner (2001) to monitor the degradation of 
hydrocarbon in the soil. This study showed a decrease in lipase activity of contaminated soil 
(Fig. 1.C and D), but the activity increased in the treated variants. This result is in agreement 
with the report of other researchers who also reported a negative correlation of lipase activity 
and reduction in TPH (Margesin and Schinner, 2001). Achieving the desired CNP ratio is an 
important consideration in bioremediation optimization (Bach et al., 2013). Microbial activity 
could be inhibited when limiting nutrients are supplied at high concentration to the 
bioremediation system, hence it is important to measure the optimum CNP ratio (Huesemann, 
1994). The CNP ratio of 100:10:1 has been reported to be the optimum ratio for oil 
bioremediation in soil (Huesemann, 1994; Wu et al., 2016). In this study, we recorded 92.95 + 
0.02% for sandy soil and 80.46 + 0.0311% for silty clay degradation with CNP ratio of 
100:2:0.2 after 36 days (Fig. 2). This indicated that there was a high reduction in TPH for 
both soils. The treatments increased soil microbial load monitored by THB, THUB, THF, 
THUF analysis (Fig. 3). Increase in the bacterial load was rapid from days 0-18 like the 
increase in fungal load from days 0-12. 

The bioremediation was successful in both sandy and silty clay soils with 93% and 
80.5% degradation of TPH respectively (Fig. 3). However, the reason for the lower TPH 
degradation rate in silty clay soil than sandy soil could be due to inefficient oxygen transfer in 
the soil. Fine grained clay with high surface area formed a sticky texture in the presence of 
water, blocking efficient oxygen transfer through the soil. Sandy soil on the other hand is 
more porous than silty clay. Higher porosity allows better oxygen transfer in the soil which is 
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essential to biodegradation of hydrocarbon. Larger pores also provide enough space for 
microbial growth. It has been demonstrated that pores smaller than 3 um are not accessible to 
bacteria (Aislabie and Deslippe, 2013). 

Another reason for low degradation of hydrocarbon in silty clay could be strong 
adsorption of the pollutant on the surface of soil particles. For the treated sandy soil, the crude 
oil residue on the 36™ day of bioremediation was 362.5 mg/kg while that of silty clay was 
418.5 mg/kg on day 36, indicating that oil remediation in sandy soil was faster compared to 
silt clay. 


Conclusion 

Crude oil pollution of soil has become recurrent in Nigeria with little or no solution 
due to over dependence on oil as the main source of revenue for the country. Therefore, 
considering the negative impact on the ecosystem, bioremediation techniques which include 
the use of fertilizers such as NPK and urea, in amounts that will not further compound the 
environmental and climate change problems, is highly recommended. 
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Abstract 

Aim of this study is to analyse the sowing qualities from six selected Robinia 
pseudoacacia L. clones from the seed orchard on the territory of State Forest Enterprise 
Parvomay (Plovdiv district). Laboratory and soil germination, energy of germination and 
sowing norms were investigated. Results show difference in the investigated indices between 
the separate clones, which probably are genetically determined. Soil germination is with about 
18% lower than the laboratory one. Clones suitable for production of black locust seedlings 
are Jaszkiseri, Appalachia and Pordim-4, whose genetic settings provide higher germination 
qualities of seeds compared to the rest of the clones. 


Keywords: black locust, clone, half-sib progeny, germination 


Introduction 

Black locust (Robinia pseudoacacia L.), which is multipurpose tree species, is 
assessed lately in many countries worldwide as one of the perspective species for production 
of biomass and energy plantations grown at higher density and short rotation (Gyuleva et al., 
2013; Gyuleva, 2014; Geyer, 2006; Rédei et al., 2010, 2011). 

The first black locust seed orchards in Bulgaria have been established in the 1980s, 
which has put the black locust management on a scientific basis (Donchev, 1989). Along with 
the development and improvement of a technology for vegetative propagation of the species 
with root cuttings (Naydenov et al., 1989; Broshtilov et al., 1998), the clonal frame for this 
species was introduced. Tsanov et al. (1992) publish the first results from testing the 
vegetative progenies of 34 selected black locust clones from 6 populations in North Bulgaria. 

From the established plantations and orchards, forest reproductive material is obtained 
today, from which seedlings with valuable qualities are produced for the needs of forestry. 

According to data from the laboratory of the Forest Seed-control Station in Sofia, 
decreasing of the vitality of black locus seed lots, originating from various regions of 
Bulgaria, is observed in the last 10 years. The literature review made for this study shows that 
there are no particular investigations of the sowing qualities of Robinia pseudoacacia L. seeds 
for the mentioned period. Similar investigations during the period from the 1950s until today 
are mainly directed to the ways of pre-sowing preparation of black locust seeds (Stefanov, 
1951; Velkov, 1968; Mirzaei et al., 2013; Basbag et al., 2010; Abdullah et al., 2019; Christin 
et al., 2019). 


Bulgarian Journal of Soil Science® 2022 Volume 7. Issue 1. www.bsss.bg 67 


Dimitrova et al., 2022 


The aim of this study is to analyse the sowing qualities of Robinia pseudoacacia L. 
seeds from six selected clones from the seed orchard on the territory of State Forest Enterprise 
Parvomay (Plovdiv district), which have shown good results according to growing indices at 
testing of their vegetative progenies in previous investigations (Tsanov et al., 1992). 


Materials and Methods 

Objects of investigation are seeds from six selected Robinia pseudoacacia L. clones, 
collected from vegetative seed orchard of State Forest Enterprise Parvomay next to Mechka 
reka nursery. The orchard was established in 2004 at a scheme 8x8 m and includes 37 clones. 
For the control, a total collection of seeds is used with provenance from vegetative seed 
orchard located in the forest nursery Lukovit (State Forest Enterprise Lesidren). 

Laboratory sowing qualities of the obtained seeds were tested in the Forest Seed- 
control Station — Sofia according to the ISTA method. 

For determining of the soil germination of clones, an experiment in the forest nursery 
Lokorsko was carried out in seed rows— in 3 repetitions at a scheme shown on figure 1. The 
seed rows includes 7 rows with a length of | m and a distance of 10 cm between the rows in 
which seeds are sown according to the scheme of each clone. Sowing was made on 15 May 
according to the norms established by the Forest Seed-control Station — Sofia. Results were 
checked every 5" day, and the last one was on 5 June. 

Pre-sowing preparation of the black locust seeds was made after the method of three 
times dipping in hot water followed by cooling. This method has been recommended as 
efficient and safe (Stefanov, 1951; Velkov, 1968; Mirzaei et al., 2013; Basbag et al., 2010; 
Abdullah et al., 2019; Christin et al., 2019). With the aim of better scarification, ice has been 
added to the cooling water. 


100 


Legend: 1 — clone Pordim 4; 2 — clone Jaszkiseri; 3 — clone Oryahovo 5; 4 — clone Appalachia; 
5 — clone Pordim 1; 6 — clone Ryahovo 7; K — control 
Fig. 1. Sowing scheme of black locust seeds in seed rows 
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Sowing was made in an universal substrate with the characteristics: medium fraction 
(0 - 40 mm), composition (in 250 1 scale: 70 % bright peat 0 - 40 mm, 30 % black peat 0 - 40 
mm, | kg NPK 14-16-18, 50 g trace elements and Tenzid — moisturising agent, pH 5,5 - 6,5. 

The sowing rate is determined according to the formula from Annex Ne7 of the 
Ordinance Ne4 from 15.02.2012 on the terms and conditions for registration of forest 
nurseries, as well as for the production of seedlings in forest nurseries — state property. 


D= 10*O*M / C*Gs*K_, where 


D — density of the sowing 

O — optimal number of ponies per 1 m ( number of annual seedlings according to 
Annex Ne7 of the Ordinance Ne4 + 20%reserve for waste) 

Gs — germination of seeds from the batch 

M — mass of 1000 seeds from the batch 

C — purity of the seeds in % 

K — correction factor for soil germination of seeds (<1, established experimentally) 


Results and Discussion 

The management of genetic resources of forest tree species and their efficient 
conservation depend on sowing qualities of seeds, which differ at various provenances of 
many species from family Fabaceae, including black locust (Rajora, Mosseler, 2001; Tozer, 
Ooi, 2014; Safdar et al., 2021; Andrea et al., 2022). 

After carrying out of analyses of sowing qualities of clone lots and the control it was 
established that the average mass of 1000 seeds varies from 17.5 to 22.7g. The laboratory 
germination is with average value 31.9% and varies from 7% to 59%. With highest 
germination are seeds from clones Jaszkiseri and Appalachia, and the first one is with 12% 
more, and the second one — with 27% lower than the control (52%). The energy of 
germination on the seventh day (Ek7) varies from 6 to 54%, and its average value is 28.7% 
(table 1). 

On the basis of physical characteristics and laboratory germination, sowing norms of 
the lots of different clones were determined; they are shown on table 1. The index values vary 
in wide range — from 0,96 to 8,25 g/m. The average value of the sowing norm is 2,7 g/m. 
Correlation between germination and seeds mass was not determined. 

Obtained results show that with highest soil germination are clones Jaszkiseri and 
Appalachia, the first one being with 14% higher and the second one — with 31% lower than 
the control (45%). The average soil germination is 26,3%. Germination of seeds begins 
between the 5th and the 10th day from the start of sowing and ends up to the 20th day. 
Percent of the germinated seeds on the 20th day varies from 5 to 52%, and its average value 
is 26,3 % (table 2). The dynamics of seeds germination, however, is different for the various 
clones and this aspect should be taken into account at the providing of reproductive material 
in connection with the final aim (Giuliani et al., 2019). 
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Table 1. Qualities of sowing materials according to clones and control in laboratory 


conditions 
Index 
Clone from Senin 
Robinia Location Piiepauee Kind of seed | Purity Abs. | Laboratory Ex7 es 
pseudoacacia | (forest enterprise) analysis mass | germination 
L. 
% G % % g/m 
1 2 3 4 5 6 7 8 9 
Pordim — 4 Vanutaliveneed Germination 
Plovdiv/ & nee after 99 19,9 31 24 
Parvomay a scarification 1,63 
Jaszkiseri Plovdiv/ Vegetative seed Germination 
Parvomay orchard after 99 19,8 59 54 
scarification 0,85 
Oryahovo — 5 Plovdiv/ Vegetative seed Germination 
Parvomay orchard after 99 17,5 21 20 
scarification 2,12 
Pordim — 1 Plovdiv/ Vegetative seed Germination 
Parvomay orchard after 99 22,7 7 6 
scarification 8,25 
Ryahovo — 7 Plovdiv/ Vegetative seed Germination 
Parvomay orchard after 99 20,7 15 13 
scarification 3,51 
Appalachia Plovdiv/ Vegetative seed Germination 
Parvomay orchard after 99 20,7 38 35 
scarification 1,39 
Control Lukovit Vegetative seed Germination 
orchard after 99 19,6 52 49 
scarification 0,96 
Average values: 99 20,1 32 28,7 2, 


It is seen from the analysis that the laboratory and soil germination of seeds is highest 


for the clones Jaszkiseri, Appalachia, Pordim-4 and the control. The applied multi-rank 
Duncan criterion identifies well-expressed and statistically significant differences in the 
average values of the laboratory germination of seeds in various clones (F=855.400, df=6, 
Sign. = 0.001), and the same one is observed in the soil germination, as well (F=890.714, 
df=6, Sign.= 0.001). 
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Table 2. Qualities of sowing materials according to clones and control in seed rows 


Index 
Date of report 
: 15 20 25 31 
Clone from ; Purity Abs. Soil ; M M M M 
Robinia Location Kind of pipe: mass | germination P ay t of ay ti “y aa 
doacacia | (forest enterprise) source _—_ acre ae ie ane nar 
pseu : Tp preparation to days; date of establishment 
: 15 May 
E15 
e .| & i ES | B10 | % | #20 
% % % 
1 2 3 4 5 6 7 8 9 10 il 
Pordim — 4 oan VeEctaye scarification 99 19,9 24 
Parvomay seed orchard 
0 19 21 24 
Jaszkiseri Plovdiv Vepciaiive scarification 99 19,8 52 
Parvomay seed orchard 
8 38 47 52 
Oryahovo — 5 Blovg Meee ave scarification 99 17,5 16 
Parvomay seed orchard 
0 9 14 16 
Pordim — 1 Bioyaty) WERetIVe scarification 99 22,7 5 
Parvomay seed orchard 
0 3 5 5 
Ryahovo — 7 Plovdiv’ See Ye scarification 99 | 20,7 11 
Parvomay seed orchard 
0 9 11 11 
. Plovdiv/ Vegetative ; : 
Appalachia Parvomay geet Orenaid scarification 99 20,7 31 
3 18 29 31 
Lesidren/ Vegetative ; : 
Control een Sead Gechail scarification 99 19,6 45 
6 36 43 45 
Average: | 99,0 | 20,1 26,3 5.7 18.9 | 24.3 | 26.3 


This presumes the existence of a correlation between the seeds provenance for both 
types of germination, which probably are genetically determined. 


In case that these clones are used for seedling production, of an interest would be also 
the testing of qualities of their half-sibs progenies with a view to obtain qualitative sowing 
material. 
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Conclusion 

The analyses of sowing qualities of seeds from selected Robinia pseudoacacia L. 
clones in the present study showed that there is a correlation between the seeds provenance 
and their laboratory and soil germination. 

Soil germination is about 18% lower than laboratory germination, which shows that 
the applied pre-sowing preparation and the selected substrate for sowing are the right ones 
and have provided optimal conditions for seeds development. 

Clones Jaszkiseri, Appalachia and Pordim-4 are suitable for the production of black 
locust seedlings, whose genetic endowments provide higher germination qualities compared 
to the rest of the clones. Robinia pseudoacacia L. is propagated both by seed and 
vegetatively. The methods of vegetative propagation may to be applied to these clones. 
Depending of the purpose for which the plantations are created, seed production of seedlings 
is much easier and cheaper method. 
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Abstract 

Historical soil science is called to analyze and evaluate realistically and critically the 
role of important dates and discoveries that determine the further development of soil science. 
Of particular importance for soil science is the first scientific and fundamental work of N. 
Pushkarov "Formation of the soil", published 112 years ago. 

With the present work an analysis of the natural and philosophical-dialectical 
interpretations of his work related to the genesis of the soil is made. Pushkarov applies the 
genetic approach to the study of soils and refutes the so-called. agro-geological (geological) 
approach, which dominated at that time in world soil science. He makes an accurate and 
objective general assessment of the factors and conditions of soil formation and proves new 
knowledge of methodical and fundamental nature, which is applied in further research in the 
country. At the same time, facts are presented about the soil science or agro-geological 
section established by him in 1911, which performs the tasks for soil-cartographic research 
formulated by him. All this gives grounds to conclude that in 1909 N. Pushkarov laid the 
foundations of soil science in Bulgaria, and in 1911 he established a scientific-administrative 
unit (scientific organization at the section level) and laid the foundations of the Institute of 
Soil Science. . 


Keywords: Nikola Pushkarov, soil science, soil science section, genetic approach, 
genesis, agro-geological bias, analysis and evaluation 


Introduction 

Two important scientific events in Bulgarian soil science, related to the name of the 
great Bulgarian natural scientist, soil scientist, public figure, patriot and revolutionary Nikola 
Pushkarov deserve to be objectively analyzed and evaluated due to the fact that they are 
extremely important for the development of soil science we have. There are scientific and 
literary sources and facts about both events, which are irrefutable proof that one event has a 
purely scientific and fundamental character, and the other - scientific and administrative. It is 
difficult to say which of them is more important, as both have their place in the development 
of Bulgarian soil science. In the order of the search for historical truth and the requirements of 
science, it is right to give priority to the first scientific event, which aims to prove the genetic 
nature of the origin of the soil. In fact, this is the first scientific work of N. Pushkarov, 
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published in the journal "Natural Science", book 8, pp. 481-489, 1909-1910, which applies 
Dokuchaev's approach and principles for scientific consideration of soil and soil formation. It 
is an indisputable fact that with this work the science of soils in our country was born and 
founded. No less important is the second event, which with the establishment at the 
suggestion of Pushkarov of a scientific-administrative organization in 1911 (soil or 
agrogeological section) at the State Agricultural Experimental Station in Sofia, laid the 
foundations of the Soil Institute in Bulgaria. With the comprehensively developed research 
program, he began to conduct the first soil surveys throughout the country. The distinction 
between these two events in role and significance is obligatory in the name of historical truth. 

What are the goals and objectives of this study? 

In connection with the above, the present study has the task to ask two important 
questions: First: to clarify the beginning of the founding of soil science in Bulgaria by 
studying the pedogenesis of Pushkarov and second: to distinguish this beginning from the 
time of the first scientific administrative unit. soil-geographical, soil-cartographic and soil- 
agrochemical studies in the country in 1911, which is incorrectly considered to be the 
birthplace of Bulgarian soil science. For this purpose, a thorough analysis and evaluation of 
the natural-scientific and philosophical-dialectical interpretations of his work, related to 
important for soil science, unexplained issues, namely: definition of soil, establishing its 
genesis, factors and conditions of soil formation, of the main and most general circular 
processes that take place in it and last but not least - upholding the genetic approach to the 
study and formation of soil. 

Historical period of writing the scientific work by Pushkarov 

Pushkarov's scientific article is still relevant today and is of lasting importance for 
both Bulgarian and world soil science. It was written when world soil science sets and 
determines the approaches, methods and directions of development of soil science. It is 
dominated by the influence of agro-cultural chemists and agro-geologists, who develop the 
agro-cultural and agro-geological bias in soil science. As a geologist-naturalist, Pushkarov 
began his first research as such and contributed to the petrography of the high Balkans 
between the peaks of Yildiz Tabia and Vezhen in 1902. As a soil scientist, however, he later 
became a successor to the soil genetic approach then unknown in us (until 1902) and founded 
in 1883 in Russia (Dokuchaev, 1883; Pushkarov, 1909). This is clearly evident in his first 
scientific work in the field of soil science and soil geological essays, which he later 
developed. 

Scientific analysis and evaluation of scientific labor 

The article, which is of fundamental importance for our science, begins with his 
famous thought "The ray gave life, and the earth - the material basis of this life." On this 
basis, he reveals some natural-genetic and his philosophical (dialectical) judgments about the 
identical nature of matter and motion, the vital motion of solar energy and earthly parts and 
the transformation of inanimate matter into living and vice versa. Describes the material 
(chemical) elements and compounds of living and non-living bodies - O, H, C, N, Fe, S, Ca, 
Mg, gases, oxides, acids and others that are the basis of soil formation of the earth's crust, 
respectively. soil and constantly enter into "new relationships" and "new groups". 

What is Pushkarov's definition of soil? In metaphorical language, he gives a very 
interesting definition of the soil, namely, "The soil is the vast laboratory where the preparation 
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of dead matter takes place in order to enter the cycle of life. Soil is a living layer of the globe. 
Like a thin, delicate diaper, it covers the vast body of the earth and gradually spiritualizes that 
body. With this definition, Pushkarov approached the doctrine of the noosphere before 
Vernadsky (2009) created it. All these interpretations so far show that he considers the soil as 
a living biological system, as part of the biosphere, which is essentially a planetary 
phenomenon of a cosmic nature. Pushkarov gives a spiritual essence to the soil. Man, as part 
of the biosphere, is the one through whom every natural body on Earth can be spiritualized. 
Man is the thinking layer of the planet Earth and his role is to cultivate or destroy soils 
(Teoharov, 2017). It is noteworthy that Pushkarov emphasizes the role of chemistry and its 
thousands of evidences for studying the cycle of elements and substances and the eternal 
transition of living to inanimate matter and vice versa. 

The role of conditions and factors in soil formation 

When describing the factors and conditions of soil formation, of course, as a geologist, 
he points out in the first place the importance of rocks in the course of soil formation 
processes, without defining them as a primary factor. For him, all factors and conditions are 
equally important in soil formation. Pushkarov believes that there is a constant exchange 
between the soil and the rock (hard crust) and the latter is the primary source of mineral 
composition and movement. He points out that the oxidation of minerals in nature is a process 
that takes place every minute. The result of this process are many metamorphic rocks and 
minerals. Therefore, the rocks form new rocks, and hence new modern soils. It takes into 
account the local nature of the rocks, which is still used as a diagnostic indicator at the generic 
taxonomic level in the classification of soils in almost all schools. According to Pushkarov, 
their local features may change under the influence of applied soils outside the location. As a 
representative example he gives the Cretaceous limestones covered with loess cover in the 
Danube plain. He considers the thick layers of rocks in the plain as former soil, which was 
carried away and deposited in the water basins gradually during the geological epochs. Later, 
our scientists confirmed the spread under the loess of the so-called. buried soils (Minkov, 
1968). Considering the question in this way, it can be concluded that much earlier than other 
scientists, in practice, he indirectly came to the theory of large and small geological 
(biological) cycle of elements and substances. In his further work in soil-cartographic 
research, Pushkarov described in detail the influence of rocks on soil formation and linked 
geological processes with soil-forming. This is proof of the extensive training he achieved in 
the early years, thanks to the fact that he had the opportunity to be guided, scientifically and 
critically evaluated by his professor, Academician Georgi Bonchev - the first Bulgarian 
petrographer and mineralogist. Along with him and other geologists, he was one of the 
founders of the Bulgarian Geological Society (Kolarova, 1959). Even today, soil science 
cannot develop without the knowledge of petrography, mineralogy and geochemistry, and 
Pushkarov applied these sciences in soil science with great competence. He considers the 
flora and fauna as important sources and factors for the formation of the soil (after the rocks). 
He called animals in the soil "organized inhabitants", noting some living organisms that are an 
integral part of it, without which it loses much of its basic character - to serve as a link 
between dead and living matter. It also emphasizes the biochemical, physical and physico- 
chemical factors of soil and nature. Pushkarov connects organisms (lower, higher), the root 
system of plants, organic residues, water and temperature with weathering (rock breaking) 
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and soil formation processes. He points out the role of water (atmospheric, rain, flowing, 
local, erosional) on weathering and soil formation, its composition and movement, 
mechanical and physical impact on minerals, rocks and organic matter. Some of his thoughts 
on its significance for the pedosphere, flora and fauna deserve attention. "The soil is the upper 
layer of the earth's crust, where the roots of plants are located and where animals find shelter 
and food. Soil that does not contain water is dead soil. It is the only solvent from which plants 
can draw mineral salts for their growth. Pushkarov describes its role in the formation of 
various geomorphological forms, and thus changes in the earth's surface, in general relief and 
topography. His example from the Balkans is remarkable, where he points out that "with its 
centuries-old work the water has cut the Balkan Mountains". He pays special attention to the 
movement of winds, emphasizing that it is the strong winds that polish and give different 
shapes to the rocks. As typical examples he points out the Belogradchik rocks, those along the 
Iskar gorge, the Stone forest near Gebedje and others. In this way he indirectly defends the 
"Aeolian hypothesis" for the formation of the most fertile fields in Europe - the Danube, 
Rhine, Yellow River and others. Pushkarov, albeit briefly, addresses the issue of 
anthropogenic activity and changes in its soil. Using folk memory, he gives a concrete 
example of the degradation of soil by man. He writes: "Old people claim that the Slivnitsa and 
Dragoman heights were covered with forests and after the destruction of the latter nothing can 
grow on the bare rock." Later, Acad. Ivan Stranski confirmed the same processes of 
deforestation and anthropogenic transformation of soils that took place after the felling of oak 
forests in the Sofia field. 

The analysis of the scientific work shows that Pushkarov strictly methodically and 
consistently applies the genetic approach in the study of genesis (soil origin). Almost all of his 
views are further embedded as an alphabetical truth in a number of works, textbooks and 
manuals. 

How does the agrogeological field of soil science study the genesis of soils at the 
same time? 

This work appears to be fundamental in the science of soils in Bulgaria. We will take 
advantage of additional facts proving its scientific and fundamental importance. To this end, 
we will make a brief comparative analysis of the approaches applied at the same time by some 
other soil schools. We will not dwell on the Dokuchaev doctrine, which, as we know, applies 
the genetic approach to the study of soils and its influence, albeit at a later stage, on world soil 
science. When Pushkarov published his work, the agro-cultural and agro-geological direction 
of studying the genesis of soils developed in Germany, which spread to a number of countries 
in Europe and America (Gerasimov, 1985). At the beginning of the twentieth century, some 
scientists in our country adhered to these approaches, without having in-depth developments 
(P. Berov - agronomist and G. Bonchev - geologist). It is important to know that agrogeology 
considers the soil as a geological formation, inert body, weathering material, which contains 
nutrients obtained from weathering of rocks. According to agrogeologists, the soil is 
considered as a weathering crust, and the process of soil formation - as a process of 
weathering (Gyurov, Kolcheva, 1969). In 1922, the Third International Conference on 
Pedology was held in Prague, attended by 14 countries, incl. and USA. It decided to detach 
soil science from the agrogeological slope and adopted the Dokuchaev approach to the study 
of soils (Stranski, 1946; Teoharov, 2016). Thus, the genetic approach to the study of 
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pedogenesis began to dominate in world soil science and factor-process, profile-evolutionary, 
horizon-profile, profile-specific (process) and evolutionary diagnostics of soils gradually 
became established; Clayden B. and JMHollis, 1984; Aubert G., 1988; Shishov et al., 2004; 
FAO, 2014). All this is proof of the relevance of Pushkarov's scientific work and its enduring 
significance. Through it he proves new knowledge of a fundamental nature, making an 
objective assessment of the main factors and conditions of soil formation that determine the 
genesis of soils and thus refutes in theory and practice the agrogeological approach in soil 
science. Thus Pushkarov, through this work, appeared as the founder of soil science in 
Bulgaria in 1909. Without making a detailed analysis, a similar opinion was expressed by 
Prof. E. Kolarova (1959). She writes: "In our opinion, this is the first article that scientifically 
examines the issues of soil and soil formation in our special literature. Moreover, this 
scientific article was written at the level of science not only at that time, but also of modern 
science - soil science. It is inexplicable why no one pays attention to this scientific assessment 
of Pushkarov's work, and all soil scientists point to 1911 as the "birth" of Bulgarian soil 
science. Another cardinal question arises from the exposition made - What was then the year 
of birth of Bulgarian soil science - 1909 or 1911? The annual report for 1910-1911 of the 
State Agricultural Experimental Station - Sofia reads as follows: "Ministry of Trade and 
Agriculture, aware of the great need for comprehensive study of soils in the country opened a 
new section at the station, very important for our agriculture agrogeological, appointing on 
17.09. the same year and a person with the appropriate preparation for its filing *(Naidenov, 
1913). This person is Nikola Pushkarov, who makes a proposal to the Minister to open a 
section on soil science. In practice, the section headed by Pushkarov is a scientific- 
administrative unit with an independent plan and tasks. In the reports on its activity he 
develops and sets for implementation the following program tasks: 1. Systematic study of 
soils by regions, for which purpose geological, biological and climatic data are collected and 
the collected samples are analyzed chemically and physically. 2. Finding out the reasons for 
the change of the properties of the soils in a given place and indicating as much as possible 
the means for elimination of the adverse phenomena. 3. Making maps for the studied areas. 
With their general formulation, these program tasks are still relevant today. These tasks give 
us reason to believe that they are strictly fundamental and applied. These are well-formulated 
tasks with a visible perspective for obtaining new knowledge and solutions for practice. They 
were adopted by the Ministry of Trade and Agriculture and their implementation began 
immediately. In 1913 the first soil-geological sketch with agrogeological map in M 1: 26000 
for the Sofia field was printed, developed independently by Pushkarov. According to Acad. 
Gerasimov (1985), the essay "is a valuable scientific monograph". In this work the opinion 
about the soil as a product of the combined influence of a number of natural factors is 
developed, which he clarifies in his work "Soil formation". Acad. Gerasimov also notes that 
this opinion is much broader and correct than previous agrogeological ideas in soil science. 


Conclusion 

"Soil formation" is the first scientific work in Bulgaria, in which Pushkarov clarifies 
the genesis and nature of soils, sets the principles of the genetic approach to their study and 
refutes the agrogeological (geological) direction in soil science. He makes an accurate and 
objective assessment of the conditions and factors of soil formation and proves new 
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knowledge of a methodical and fundamental nature, which still has its current relevance. With 
the developed program tasks and the soil science section created on his proposal, he applies 
this approach and the new knowledge in science and practice. All this is proof that Pushkarov 
founded soil science in 1909, and through the scientific and administrative organization he 
created at the section level in 1911, he laid the foundations of the Institute of Soil Science. 
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cm) which is the maximum size for the illustrations. Magnification should be shown by scale 
bars. Files in XLS, TIF, JPG and WME formats are accepted. After the paper is accepted the 
editors will ask for figures in JPG or TIF format with 300 Dpi. Colour illustrations are 
accepted without fees for electronic publications. The editors reserve the right to decide what 
figures are justified and can be published in RGB colour (other figures will be removed or 
published in black and white). 

Example: Fig. 1. Soil map of the study area. 
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Psychological Association) bibliographic style (7th edition — 2020). APA Referencing Basics: 
In-Text Citation. 

Citation 
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References should be arranged first alphabetically and then further sorted 
chronologically if necessary. More than one reference from the same author(s) in the same 
year must be identified by the letters ‘a’, ‘b’, ‘c’, etc., placed after the year of publication. 
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